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. Motivated by the problem of the dynamics of point-particles in high post- 

. Newtonian (e.g. 3PN) approximations of general relativity, we consider a cer- 

O , tain class of functions which are smooth except at some isolated points around 

I , which they admit a power- like singular expansion. We review the concepts 

of (i) Hadamard "partie finie" of such functions at the location of singular 
points, (ii) the partie finie of their divergent integral. We present and inves- 
^ , tigate different expressions, useful in applications, for the latter partie finie. 

To each singular function, we associate a partie-finie (Pf) pseudo-function. 
The multiplication of pseudo-functions is defined by the ordinary (pointwise) 
product. We construct a delta-pseudo- function on the class of singular func- 
tions, which reduces to the usual notion of Dirac distribution when applied 
on smooth functions with compact support. We introduce and analyse a new 
derivative operator acting on pseudo- functions, and generalizing, in this con- 
text, the Schwartz distributional derivative. This operator is uniquely defined 
up to an arbitrary numerical constant. Time derivatives and partial deriva- 
tives with respect to the singular points are also investigated. In the course of 
the paper, all the formulas needed in the application to the physical problem 
are derived. 
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I. INTRODUCTION 



The Hadamard regularization 1T|,0], based on the concept of finite part ("partie finie") 
of a singular function or a divergent integral, plays an important role in several branches of 
Mathematical Physics (see |^|-§[ for reviews). Typically one deals with functions admitting 
some non-integrable singularities on a discrete set of isolated points located at finite distances 
from the origin. The regularization consists of assigning by definition a value for the function 
at the location of one of the singular points, and for the (generally divergent) integral of 
that function. The definition may not be fully deterministic, as the Hadamard partie finie 
depends in general on some arbitrary constants. The Hadamard regularization is one among 
several other possible regularizations 

A motivation for investigating the properties of a regularization comes from the physical 
problem of the gravitational interaction of compact bodies in general relativity. As it is 
hopeless to find a sufficiently general exact solution of this problem, we resort to successive 
post-Newtonian approximations (limit c — * +00). Within the post-Newtonian framework, 
it makes sense to model compact objects like black holes by point-like particles. This is 
possible at the price of introducing a regularization, in order to cure the divergencies due to 
the infinite self-field of the point-masses. However, general relativity is a non-linear theory 
and, if we want to go to high post-Newtonian approximations, involving high non-linear 
terms, the process of regularization must be carefully defined. In particular, it turns out 
that, from the third-post-Newtonian approximation (3PN or 1/c^), the problem becomes 
complicated enough that a rather sophisticated version of the Hadamard regularization, 
including a theory of generalized functions, is required. By contrast, a cruder form of 
the Hadamard regularization, using merely the concept of partie finie of singular functions 
l^r^l, is sufficient to treat the problem up to the 2PN order. Furthermore, we know 
that the answer provided by the Hadamard regularization up to the 2PN order is correct, 
in the sense that the field of the two bodies matches the inner field generated by two 
black holes |]14|, and the result for the equations of motion can be recovered without need 



of any regularization from computations valid for extended non-singular objects |T5,16 



Conforted by these observations we systematically investigate in this paper the Hadamard 
regularization as well as a theory of associated generalized functions, in a form which can 
be directly applied to the study of the dynamics of two point-like particles at the 3PN order 
]T7| . (We therefore restrict our attention to two singular points; however most of the results 
of the paper can be generalized to any number of points.) Notice that this problem enjoys 
a direct relevance to the future gravitational-wave experiments LIGO and VIRGO, which 
should be able to detect the radiation from black-hole and/or neutron-star binaries which a 



precision compatible with the 3PN approximation |T8 



Consider the class JF of functions on that are smooth except at two isolated singulari- 
ties 1 and 2, around which they admit some power-like singular expansions. The Hadamard 
partie finie (F)i of F G at the location of singularity 1, as reviewed in Section II, is 
defined by the average over spatial directions of the finite-part coefficient in the expansion 
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of F around 1. On the other hand, the Hadamard partie finie Pf J d^xF of the divergent 
integral of F, we will review in Section III, is obtained from the removal to the integral 
of the divergent part arising when two regularizing volumes surrounding the singularities 
shrink to zero. Both concepts of partie finie are closely related. Notably, the partie-finie 
integral of a gradient is equal to the sum of the parties finies (in the former sense) of the 
surface integrals surrounding the singularities, in the limit of vanishing areas. In Section IV 
we investivage several alternative expressions of the Hadamard partie finie of integrals, some 
of them based on a finite part defined by means of an analytic continuation process (see 
for a relation between partie finie and analytic continuation). In our terminology, we adopt 
the name "partie finie" for the specific definitions due to Hadamard, and speak of a "finite 
part" when referring to other definitions, based for instance on analytic continuation. In 
Section V we focus to the case (important in applications) of the partie finie of a Poisson 
integral of F G JF. 

To any F G JF, we associate in Section VI a generalized function, or partie-finie "pseudo- 
function" PfF, which is a linear form on JF defined for any G G by the duality bracket 
< PfF, G >= Pf j d^:KFG. When restricted to the set T> of smooth functions with compact 
support the pseudo-function PfF is a distribution in the sense of Schwartz (see also 
||T9|^I1 for more details about generalized functions and distributions), i.e. a linear form 
which is continuous with respect to the Schwartz topology. [However, we do not attempt 
here to introduce a topology on JF; we simply define the set of algebraic and differential rules, 
needed in applications, that are satisfied by the pseudo-functions on JF.] The product of 
pseudo-functions coincides with the ordinary ( "pointwise" ) product used in Physics, namely 
PfF . PfG = Pf(FG). An important particular case is the pseudo-function Pf^i obtained 
(in Section VI) from the pseudo-function associated with the Riesz delta-function ||22|, and 
that satisfies VG G JF, < Pi6i,G >= (G)i. The "Dirac pseudo-function" Pf5i plays in 
the present context the same role as plays the Dirac measure in distribution theory. We 
introduce also more complicated objects such as Pf(F5i). In Sections VII and VIII we 
show how to construct a derivative operator on JF, generalizing for this class of function 
the standard distributional derivative operator on V and satisfying basically the so-called 
rule of integration by parts, namely VF, G G JF, < (9i(PfF),G >= - < 9i(PfG),F >. In 
addition we require that the derivative reduces to the "ordinary" derivative for functions that 
are bounded in a neighbourhood of the singular points, and that the rule of commutation 
of derivatives holds. We find that this derivative operator is uniquely defined modulo a 
dependence on an arbitrary numerical constant (see Theorem 4 in Section VIII). It represents 
a natural notion of derivative within the context of Hadamard regularization of the functions 
in JF. However, it does not satisfy in general the Leibniz rule for the derivative of a product 



(in agreement with a theorem of Schwartz ||2^). See Colombeau ||2J] for a multiplication of 
distributions and associated distributional derivative satisfying the Leibniz rule. Further, 
we obtain the rules obeyed by the new derivative operator when acting on pseudo-functions 
such as Pf(F5i) in Section VII, and we investigate the associated Laplacian operator in 
Section VIII. Finally, in Section IX, we consider the case of partial derivatives with respect 
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to the singular points 1 and 2, as well as the time derivative when both singular points 
depend on time (i.e. represent the trajectories of real particules). Within this approach, the 
latter distributional derivative constitutes an important tool when studying the problem of 
the gravitational dynamics of point-particles at the 3PN order . 



Notation. N, Z, M and C are the usual sets of non-negative integers, integers, real numbers 
and complex numbers; M"*"* is the set of strictly positive real numbers s > 0; M'^ is the usual 
three-dimensional space endowed with the Euclidean norm |x| = {xl + + Xg)^/^; C^{il) is 
the set of p-times continuously different iable functions on the open set Q {p < +oo); Ll^^{fl) 
is the set of locally integrable functions on f2; the o and O symbols for remainders have 
their standard meaning; distances between the field point x and the source points yi and 
y2 are denoted by ri = |x — yi| and r2 = |x — y2|; unit directions are ni = (x — yi)/ri 
and n2 = (x — y2)/r2', dQi and dQ2 are the solid angle elements associated with rii and n2; 
^12 = |yi — y2|; 'Si(s) and -62(5) denote the closed spherical balls of radius s centered on yi 
and y2; (9, = d/dx\ idi = d/dyl, 2<9j = d/dy\] L = iii2 ■ ■ -i; is a multi-index with length 
/; = n'^ ■ ■ -n^ and = di^- ■ ■ d^; the symmetric-trace-free (STF) projection is denoted 
by = STF(?T,f'); (ij) = ^^^^ and [ij] = ^^^^^ 1 2 means the same expression but 
corresponding to the point 2; for clearer reading, we use left-side labels 1 and 2 when the 
quantity appears within the text, like for the partial derivatives idi and 2di or the coefficients 
ifa and 2/b, and labels placed underneath the quantity when it appears in an equation; iff 
means if and only if. 



II. HADAMARD PARTIE FINIE 

A. A class of singular functions 

All over this paper we consider the class of functions of a "field" point x G that are 
singular at the location of two "source" points yi and y2 around which they admit some 
singular expansions. 

Definition 1 A real function -F(x) on is said to belong to the class of functions T iff 
(i) F is smooth on M.^ deprived from yi and y2, i.e. F G C°°(]R'^ — {yi,y2})- 
(a) There exists an ordered family of indices (ai)igN with G M, and a family of coeffi- 
cients ^fa^, such that 

iff 

ViV G N , F(x) = J]rr /a.(ni)+ Rn{^) . (2.1) 

i=0 

Here ri = |x — yi| and ni = (x — yi)/ri; i^ satisfies Oq < o-i < ■ ■ ■ < < N < Oij^+i; and 
the "remainder" is 

i^^(x) = o(rf ) when n ^ . (2.2) 
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(Hi) Idem with indices (&i)iGN; coefficients ^fb^, remainder 2Rn, ri ^ r2 and rii ^ n2. 

In addition to Definition 1, we always assume that the functions F & J-' decrease sufficiently 
fast at infinity (when |x| oo) so that all integrals we meet are convergent at infinity. 
Thus, when discussing the integral / d^^KF, we suppose implicitly that F = o(|x|^^) at 
infinity [or sometimes F = 0{\x.\~^~^) where e > 0], so that the possible divergencies come 
only from the bounds at the singular points yi,2- Similarly, when considering the integral 
f d^yiFG, we suppose FG = o(|x|~'^), but for instance we allow F to blow up at infinity, say 
F = 0(|x|), if we know that G decreases rapidly, e.g. G = o(|x|~^); in the case of / d^xc^jF, 
we generally assume F = o(|x|^^). [Clearly, from Definition 1 the ordinary product FG of 
two functions of JF is again a function of JF; and similarly the ordinary gradient diF G J-".] 

An important assumption in Definition 1 is that the powers of ri in the expansion of 
F when Vi —>■ (and similarly when r2 0) are bounded from below, i.e. Oq < ctj where 
the most "divergent" power of ri, which clearly depends on F, is oq = ao{F). Thus the 
part of the expansion which diverges when ri ^ is composed of a finite number of terms. 
Notice also that we have excluded in Definition 1 the possible appearance of logarithms 
of Ti (or in the expansion of F. See Sellier [^ for a more general study in the case 
where some arbitrary powers of logarithms are present. We will discuss the occurence of 
logarithms in Section V, when dealing with the Poisson integral of F. At last, we point out 
that the coefficients ifa (and similarly 2 ft) do not depend only on rii, but also they do on 
the source points yi and y2, so that in principle we should write i/a(ni; yi, y2); however, 
for simplicity's sake we omit writing the dependence on the source points. The coefficients 
could also depend on other variables such as the velocities vi and V2 of the source points, 
but the velocities do not participate to the process of regularization and can be ignored 
for the moment (we will return to this question in Section IX when considering the time 
dependence of F). 

Once the class has been defined, we shall often write in this paper the expansions of 
F when ri_2 — in the simplified forms 

F(x)= /„(ni) + o(rf ) when n , (2.3a) 

ao<a<N 

^(x)= YI 4fbM + o{r^) whenra^O, (2.3b) 

bo<b<N 

by which we really mean the expansions in Definition 1, i.e. in particular where the indices 
a G (aj)ieN and b G (&i)jgN, and are a priori real. However, most of the time (in applications), 
it is sufficient to assume that the powers of ri 2 are relative integers a, 6 G Z. We can then 
write the expansion ri ^ in the form 

fco N 

P = llzm + E ^1 + ' (2-4) 

k=0 1 k=0 
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where fco = — 1 — ao- In the following we shall sometimes derive the results in the simpler case 
where the powers G Z, being always undertood that the generalization to the case of real 
powers is straightforward. Finally, it is worth noting that the assumption (i) in Definition 
1, that F is C°° outside {yi,y2}, can often be relaxed to allow some functions to have 
integrable singularities. An example is the function x — > l/|x — x'| encountered in Section 
V, depending on a fixed "spectator" point x' distinct from yi and y2. To treat such objects, 
we introduce a larger class of functions, J-'ioc- 

Definition 2 -F(x) is said to belong to the class of functions T\oc iff 

(i') F is locally integrable on deprived from yi and y2, i.e. F G Lj'Q^(M^ — {yi,y2})- 
(a)- (Hi) in Definition 1 hold. 

For simplicity, in the following, we shall derive most of the results for functions belonging 
to the class (even if the generalization to J^\oc is trivial); J-'ioc will be employed only 
occasionally. 

B. Partie finie of a singular function 

The first notion of Hadamard partie finie is that of a singular function at the very location 
of one of its singular points. 

Definition 3 Given F & J-' we define the Hadamard partie finie of F at the point yi to be 

(F), = I ^ foM , (2.5) 

where dfli = dfl{ni) denotes the solid angle element of origin yi and direction ni. 

In words, the partie finie of F at point 1 is defined by the angular average, with respect to 
the unit direction ni, of the coefficient of the zeroth power of ri in the expansion of F near 
1 (and similarly for the point 2). There is a non zero partie finie only if the family of indices 
(aj)jgN in Definition 1 contains the value 0, i.e. 3zo such that = 0. The latter definition 
applied to the product FG of two functions in JF yields 

ao(F)<a<-ao(G) 

where if a and iQa are the coefficients in the expansions of F and G when ri — (the 
summation over a is always finite). From ( |2.6| ) it is clear that the Hadamard partie finie is 
not "distributive" with respect to the multiplication, in the sense that 

{FG)i^{FUG)i in general. (2.7) 

The partie finie picks up the angular average of i/o(ni), namely the scalar or / = piece 
in the spherical-harmonics expansion [Yim), or, equivalently, in the expansion on the basis 
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of symmetric and trace-free (STF) products of unit vectors ni = {n\). For any / G N, we 
denote by L = iii2 ■ ■ - ii a. multi-index composed of / indices, and similarly L—1 = 11X2 ■ ■ ■ ii-i, 
P = iij2 ■ ■ ' jp- In general we do not need to specify the carrier index i or j, so a tensor with / 
upper indices is denoted T^, and for instance the scalar formed by contraction with another 
tensor of the same type is written as S* = T^U^ = T*^'"*'f/*i '*', where we omit writing 
the / summations over the I indices ik = 1, 2, 3. We denote a product of I components of the 
unit vector n\ by = n^^ ■ ■ -n-l, and the STF projection of that product by = STF(nf'): 
e.g., hj' = n\n{ — '^6^\ n^^ = n\n{n^ — ^{n\6^'^ + n{6^^ + n^6^^). More generally, we denote by 
the STF projection of T^; that is, is symmetric, and satisfies 5j j T*'-i*'^~^ = (see 



| 25| and the appendix A of for a compendium of formulas using the STF formalism). 



The coefficients if a of the expansion of F admit the STF decomposition 

+00 

/,(nO = 5^ f/^ , (2.8) 

1=0 

where the i/i^'s are constant STF tensors, given by the inverse formula: 



,1 (2(+l)!! rddi 



In STF notation, the Hadamard partie finie of F at 1 reads simply 

{F)i = fo , (2.10) 
where if a denotes the first term in the expansion (p.8|). 



Lemma 1 The partie finie at 1 of the gradient diF ( as defined outside the singularities ) of 
any function F ^ T satisfies 

(diF)i = zC^F\ . (2.11) 
\ri /I 

This Lemma is particularly useful as it permits replacing systematically the differential 
operator di by the algebraic one 3^ when working under the partie-finie sign (..)i. 

Proof. The expansion when ri — of the gradient is readily obtained from the expansion 
of F itself as 

d,F = J2 [« < fa + difa] , (2.12) 

a 

(with over-simplified notation for the sum), where the operator dl is defined as ridi when 
applied on a function of the sole unit vector iii. Hence, explicitly, dl = {5^^ — ni)-^- This 

operator is evidently transverse to rii : n\dl = 0, and we get, from the decomposition ( p.SD , 
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1=0 

Thus, by averaging over angles, 



+00 

dl fa = Yl ^(^1"' f"^"' - ^ (a) ■ (2-13) 



''''dlU = lf: = 2['4^n\f.. (2.14) 



47r ' I 3 r 471 ' 1 



We readily deduce that the partie finie of the gradient (|2.12| ) is given by 

= 3l^n\ A = p (QED). (2.15) 

As an example of application of Lemma 1, we can write, using an operation by parts, 
(rf(?iF)i = [dii^rfF) - di{rl)F]i = [3n\rlF - di{rf)F]i, from which it follows that 

{rldiF), = 0. (2.16) 

Another consequence of Lemma 1, resulting from two operations by parts, is (rf AF)i = 
[3n\ridiF — di{rl)diF]i = {n\ridiF)i = [3F — di{n\ri)F]i (where the Laplacian A = 9,9,), 
hence the identity 

(r^AF)! = . (2.17) 

By the same method we obtain also 

{d^,F)^ = ( '"^^ , f) = 2 ( f,'^ + 5^^ f,) , (2.18) 
\ ri /II 1 

the right-hand side of the last equality being expressed in terms of the STF tensors 
parametrizing (|2.8| ). Tracing out the previous formula, we find 

(AF)i=(-|f) =6/2. (2.19) 

Vrf /I 1 

Finally, let us quote the general formula for the partie finie of the Zth derivative d^F = 
di^ ■ --di^F: 

(9lF)i = UY, ^^'"^ //""'"^^ • (2.20) 

fc=0 ^ 

Here, [|] denotes the integer part of |, S"^^ is the product of Kronecker symbols 
§iii2§i3i4 . . .§i2k-ii2k^ q^t^^ ^/^^"^^ = i/;*^*"^^'"*' ; the parenthesis around the indices denote the 
symmetrization. One may define the "regular" part of the function F near the singularity 
1 as the formal Taylor expansion when ri — obtained using ( |2.2(]|) . Thus, 

+00 ^ +00 L2J 

^1" - E lr[n'^i9,F)^ = J2r[i: n{-'^ f/''^ • (2.21) 

1=0 ' 1=0 k=0 
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III. PARTIE-FINIE INTEGRALS 

A. The partie finie of a divergent integral 

The second notion of Hadamard partie finie is that of the integral J (i^xF(x), where 
F E J-". This integral is generally divergent because of the presence of the singular points 
yi and y2 (recall that we always assume that the function decreases sufficiently rapidly at 
infinity so that we never have any divergency coming from the integration bound |x| +00). 
Consider first the domain deprived from two spherical balls Bi{s) and B2{s) of radius s, 
centered on the two singularities yi, Bi{s) = {x; ri < s} and 82(3) = {x; r2 < s}. We 
assume that s is small enough, i.e. s < ^ where ri2 = |yi — y2|, so that the two balls do 
not intersect. For s > the integral over this domain, say J(s) = /iR3\g^(s)uB2{s) '^^''^ 
well-defined and generally tends to infinity when s — 0. Thanks to the expansions (assumed 
in Definition 1) of F near the singularities, we easily compute the part of I{s) that blows up 
when s — s> 0; we find that this divergent part is given, near each singularity, by a finite sum 
of strictly negative powers of s (a polynomial of 1/s in general) plus a term involving the 
logarithm of s. By subtracting from I{s) the corresponding divergent part, we get a term 
that possesses a finite limit when s —>■ 0; the Hadamard partie finie H is defined as this 
limit. Associated with the logarithm of s, there arises an ambiguity which can be viewed as 
the freedom in the re-definition of the unit system we employ to measure the length s. In 
fact it is convenient to introduce two constant length scales si and S2, one per singularity, 
in order to a-dimensionalize the logarithms as In(^) and In(^). 

Definition 4 For any F E J-' integrable in a neighbourhood of |x| = +00, we define the 
Hadamard partie finie of the divergent integral J d^x F as 

Pf,„,, [d'^F = hm { / rf^x F + V ^ / dn^ + In ( -] [ dU, /_3 

+ 1^2}, (3.1) 

where 1^2 means the same previous two terms but concerning the singularity 2. 

This notion of partie finie can be extended to functions which are locally integrable outside 
the singularities, i.e. F G T\oc (see Definition 2). In ( ^.1| ) the divergent terms are composed 
of a sum over a such that a + 3 < as well as a logarithmic term, by which we really mean, 
using the more detailed notation of Definition 1, 

E ^ / . h + ) / + 1 " 2 ■ 

where ii is such that < a\ < ■ ■ ■ < ctjj-i < —3 < (the sum is always finite); we 
have introduced a Kronecker symbol 5-3,0^^ to recall that the logarithm is present only if 
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the family of indices (ai)jgN contains the integer —3 (i.e. = —3). The divergent terms in 
can also be expressed by means of the partie finie defined by (|2.5|). Indeed, they read 



+ 1^2 



a+3<0 

[coming back to the less detailed notation of (|3.1|) j. 

The partie-finie integral (|3.1|) depends intrinsically on the two arbitrary constants si 
and S2 introduced above. There is another way to interpret these constants besides the 
necessity to take into account the dimension of s, which is discussed by Sellier in 0. With 
this point of view we initially define the partie finie using two arbitrarily shaped volumes Vi 
and V2 instead of the two spherical balls Bi and 82- Consider for instance the two volumes 
Vi = {x; ri < spi(ni)} and V2 = {x; r2 < sp2(n2)}, where s G M"*"* measures the size 
of the volumes and the two functions pi and p2 describe their shape (the balls Bi and B2 
corresponding simply to pi and p2 = 1). Here, we assume for simplicity that the volumes 
remain isometric to themselves when s varies. Then, the partie finie is defined as the limit 
of the integral over \ Vi U V2 to which we subtract the corresponding divergent terms 
when s — > 0, without adding any normalizing constant to the logarithms. In this way, we 
find that the alternative definition is equivalent to our definition (|3.1|) provided that si and 
S2 are related to the shapes of the regularizing volumes Vi and V2 through the formula 

Insi j dn, /_3 = J Alnpi , (3.2) 

(and similarly for S2). The arbitrariness on the two original regularizing volumes is therefore 
encoded into the two (and only two) constants si and S2. A closely related way to interpret 
them is linked to the necessity to allow the change of the integration variable x in the integral 
J (i^xF. Such an operation modifies the size and shape of the regularizing volumes, thus the 
balls Bi and B2 are in general transformed into some new volumes Vi and V2; so, according 
to the previous argument, the freedom of choosing the integration variable reflects out in 
the freedom of choosing two arbitrary constants Si and S2. (In this paper we shall assume 
that Si and S2 are fixed once and for all.) 

An alternative expression of the Hadamard partie finie is often useful because it does 
not involve the limit s — 0, but is written with the help of a finite parameter s' G M"*"*. 
Consider some s' such that < s < s', and next, split the integral over R^\Si(s) Ufi2(s) into 
the sum of the integral over M.^ \ Bi{s') U -62(5') and the two integrals over the ring-shaped 
domains Bi{s') \ Bi{s) and Bi{s') \Bi{s). If s < s' <^ 1 we can substitute respectively into 
the ring-shaped integrals the expansions of F when ri — > and r2 [see ( p.3|) ]. The 
terms that are divergent in s cancel out, so we can apply the limit s — > (with fixed s'). 
This yields the following expression for the partie finie: VA^ G N, 

) / dn, /_3 



pf,„,, / d'^F= [ d'^F+ V ^ / dn,f, + \n(-] [ dn,f^ 



a + 3#0 



+ 1^2 + o(s''^+') , (3.3) 
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which is vahd for an arbitrary fixed s'. Of course, up to any given finite order N the second 
member of ( p.3| ) depends on s', but in the formal hmit +00, this dependence disappears 
and, in fine, the partie finie is independent of s'. 

B. Partie-finie integral of a gradient 

A fundamental feature of the Hadamard partie finie of a divergent integral is that the 
integral of a gradient diF is a priori not zero, since the surface integrals surrounding the 
two singularities become infinite when the surface areas shrink to zero, and may possess a 
finite part. 

Theorem 1 For any F E J-" the partie finie of the gradient of F is given by 

Pf j rf^x d,F = -4:71 in\rlF)i + 1^2, (3.4) 
where the singular value at point 1 is defined by (|2.5|) . 



In the case of a regular function, the result is always zero from the simple fact that the 
surface areas tend to zero — c/the factor r\ in the right side of (|3.4| ). However, for F G JF, 
the factor r\ is in general compensated by a divergent term in the expansion of F, possibly 
producing a finite contribution. 



Proof. We apply ( |3.1D to the case of the gradient d^F, using the expansion of diF when 
ri — s> as given by ( p.l2| ). The expression of the divergent terms is simplified with the help 
of the identity ( |2.14| ), which shows notably that the logarithms and associated constants $1^2 
disappear. This leads to 



lim 



R3\Bi{s)uB2(s) 



d"^^ d,F + I dn^nl + 1 ^ 2| . (3.5) 



a+2<0 



Next, the first term inside the braces is transformed via the Gauss theorem into two surface 
integrals at ri = s and r2 = s, where we can replace F by the corresponding expansions 
around yi and y2 respectively. We get 

lim s"^' I < fa + / < f^} = ~ I "^^1 h 

a a+2<0 

(and similarly when 1^2); QED. 

From Theorem 1 it results that the correct formula for "integrating by parts" under the 
sign Pf is 

Pf j d^x FdiG = -Pf j d^x GdiF - Air {n\rlFG)i - An {v^^rlFG)2 . (3.6) 
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Note also that the partie-finie integrals of a double derivative as well as a Laplacian are 
given by 

Pf J rf^x dijF = 4n (^ri{6'^ - 2<)f)^ + 1^2, (3.7a) 

Pf J £x AF = 47r (riF)i + 1 2 . (3.7b) 

C. Parties finies and the Riesz delta-function 

The Riesz delta-function plays an important role in the context of Hadamard parties 
finies. It is defined for any e G M.^* by e5(x) = ^^i^|x|^~^; when £ — >■ 0, it tends, in the 
usual sense of distribution theory, towards the Dirac measure in three dimensions — i.e. 
limg^oe^ = as can be seen from the easily checked property that A(|x|^^^) = —An e5(x). 
The point for our purpose is that when defined with respect to one of the singularities, the 
Riesz delta-function belongs to JF. Thus, let us set, We G M+*, 

,5i(x) = ,5(x - yi) = G ^ (3.8) 

An 

(and idem for 2). Now we can apply to e5i(x) the previous definitions for parties finies. In 
particular, from Definition 3, we see that has no partie finie at 1 when e is small enough: 
{e^i)i = 0. From Definition 4: 

Lemma 2 For any F E T , we have 

lim Pf / rf^x = (F)i , (3.9) 
where the value of F at point 1 is given by the prescription ( |2.5| ) . 

Proof. For e > we evaluate the finite part of the integral for the product i^SiF G T using 



the specific form (|3.3| ) of the partie finie defined in terms of a given finite s' . The expansions 



of e^iF when ri^2 tend to zero are readily determined to be 

e5,F = iiizfl Vrr^-Va(ni) for n ^ , (3.10a) 

47r ^ — ' 1 

a 

= £ilz£) V ^ 9Lrf2-'Vr^+'n^/,(n2) for ^ . (3.10b) 

An ^-^ II 1 ^-^ 2 

/>0 b 

In the second equation we used the Taylor expansion rf~^ = J2i>o ■^7!^^2^2 i'9i'"i2^^ when 
r2 — 0, with notation n2 = n'2'''^2 idi = idi^ ■ ■ ■ id^. Hence, we can write the 
partie-finie integral in the form (VA^ G N; with fixed s' such that < s' < 1) 
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i>0 6+i+3<JV 

Here, we have discarded the term with In (^j^^ by choosing e > to be so small that all 
denominators a + e differ from zero. Since tends towards the Dirac measure when e 0, 
the integral over M'^ \ Bi{s') U -62(5') goes to zero. Because of the factor e present in the 
numerators, so do the other terms when e —>■ 0, except for those whose denominators involve 
a compensating e. Now, the only term having the required property corresponds to a = 
in the previous expression. Therefore, taking the limit e — (with fixed s'), we get 

lim Pf J rf^x ,6,F = J ^/o(ni) + o(.'^) , 

and this being true for any A^, we conclude 

lim Pf J d^x ,5,F = J ^/o(ni) = (F), (QED). 

As we can infer from Lemma 2, the Riesz delta-function should constitute in the limit 
e — >■ an appropriate extension of the notion of Dirac distribution to the framework of parties 
finies of singular functions in JF. The precise definition of a "partie-finie Dirac function" 
necessitates the introduction of the space of linear forms on and will be investigated in 
Section VI (see Definition 7). 

IV. ALTERNATIVE FORMS OF THE PARTIE FINIE 
A. Partie finie based on analytic continuation 

Practically speaking, the Hadamard partie-finie integral in the form given by ( ^.1| ) is 
rather difficult to evaluate, because it involves an integration over the complicated volume 

\ Bi{s) U 82(3). Fortunately, there exist several alternative expressions of the Hadamard 
partie finie, which are much better suited for practical computations. The first one is based 
on a double analytic continuation, with two complex parameters a, G C, of the integral 



.^x(-il l-i) F, (4.1) 



where the constants Si and S2 are the same as those introduced within the definition ( ^.ip . 
The point for our purpose is that the integral ([4.1|) does range over the complete set R^. First 
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of all, we propose to check that Ia,i3 is defined by analytic continuation in a neighbourhood 
of the origin a = = /3 in C^, except at the origin itself where it generically admits a simple 
pole in a or /3 or both. We start by splitting J^,/? into three contribution: ila,f3 extending over 
the ball Bi{s) of radius s surrounding 1, 2-^a,/3 extending over the ball 82(3) surrounding 2, 
and 3la,i3 extending over the rest M^\i3i(s) US2(s). The integral ila,i3 is initially convergent 
for 3fJ(a) > — ao — 3 and any (3, where ao is the most singular power of ri in the expansion of 
F near yi; similarly, 2-^a,/3 exists only if 3fJ(/3) > —60 — 3 and any a {bo is the analogous to oq 
that relates to y2), and 3/o,/3 exists if ^{a + [3) < e, where e > is such that F = 0(|x|^^~^) 
when |x| +00. As the third contribution 3/0,/? is clearly defined in a neighbourhood of 
the origin, including the origin itself, we consider simply the part i/o,/3 (the same reasoning 
applies to 2-^0,/?) • Within the integrand, we replace the product r2-F by its expansion in the 
neighbourhood of yi (using a Taylor expansion for rg), and find that the dependence on 
/3 occurs through some everywhere well-defined quantity, namely i^^rfg- After performing 
the angular integration over dQi, we obtain a remaining radial integral consisting of a sum 
of terms of the type drir""'"'^"'''"^^ = (^q, _|_ _|_ / _|_ 3)^ ^^j^at clearly admit a unique 

analytic continuation on C\Z; hence our statement (a simple pole at the origin arises when 
a = -/-3). 

Theorem 2 For any function F E J-" that is summahle at infinity, the Hadamard partie 
finie of the integral is given by 

where FP a^o means taking the finite parts in the Laurent expansions when a — > and (3 —>■ 
successively. 

The proof of Theorem 2 is relegated to Appendix A. Notice our convention regarding the 
notation: while "Pf" always stands for the Partie finie of an integral in the specific sense of 
Hadamard 0] , we refer to "FP" as the Finite Part or zeroth-order coefficient in the Laurent 
expansion with respect to some complex parameter {a, /5 G C, or i? G C as in the next 
subsection). We see from Theorem 2 that the partie finie Pf can be viewed as a finite part 
FP and vice versa. The link between analytic continuation and Hadamard partie finie is 
pointed out by Schwartz 0]. More precisely. Theorem 2 says how to calculate the Hadamard 
partie finie; the procedure consists of: (i) performing the Laurent expansion of Ia,i3 when 
a — > while jS remains a fixed ("spectator") non-zero complex number, i.e. 

+00 

P=Pmin 

where p gZ and where the coefficients /(p),/3 depend on f3; (ii) achieving the Laurent expan- 
sion of the zeroth-a-power coefficient /(o),/3 when /? — >• 0, i.e. 

+00 

^(0),/3 = Ao,g) , 

?=<?min 
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to finally arrive at the zerotli-/3-power coefficient /(o,o)- Indeed, we find that the same result 
can be obtained by proceeding the other way around, first expanding around (3 = with a 
fixed a, then expanding the coefficients Ia,{o) near a = 0. Thus, 



fp|fp/„,^ 

/3-»0 Iq!-+0 



'(0,0) 



:Fp|fp/„,4 • 



(4.3) 



We emphasize that the definition ( |3.1| ) of the partie finie yields unambiguously the result 
1(0,0), which corresponds to taking independently the two limits a — >• and P ^ (the 
limiting process does not allow for instance to keep a = (3). The final value /(o,o) is the same 
as the one given by the regularization adopted by Jaranowski and Schafer []T2| (see their 
appendix B.2). 

In practice the expression ( |4.2|) is used in connection with the Riesz formula |2^, valid 
for any 7, 5 G C except at some isolated poles. 



3/2 



r (i±5) F (^3) r (- 



7+^+3 



rH)r(-|)r(^) 



■y+5+3 
12 



(4.4) 



with = |yi — y2|; here, F denotes the Eulerian function. According to Theorem 2, the 
formula ( [4. 4] ) permits computing the partie finie of any integral of a product between powers 
of ri and r2. Consider the (not so trivial) case of the integral of rjf'^r^^, which is divergent 
at both points 1 and 2. From the Riesz formula, with 7 = a — 3 and 5 = /? — 3, we have 



TT' 



3/2 _ 



r(|)r(f)r( 



a+/3-3 



) 



r 



a+l3-3 
12 



r(-^)F(-^)F(^) sfs 



We compute the Laurent expansion when a ^ with fixed /3 G C and obtain a simple pole 
in a followed by a /5-dependent finite part given by 



'(0),/3 



TT 



3/2 



F(l) r 



(3-3 
12 



+ vl/(l)-vl;(l + f)+vl.(|)_vl;(|_|)+21n 



ri2 

Sl 



with "^(z) = ^ In T(z). This finite part itself includes a simple pole in P, and then we obtain 
the corresponding finite part when /3 — as 



'(0,0) 



n'/^ F(l) 

r'u r (I) 

At last. Theorem 2 tells us that 



21n^ 



Pf. 



4lT 



3 ^3 



r\ r. 



1 ' 2 



12 L 



n_2 

Sl 



ri2 



+ 21n(^ 

S2 



Sl 



In ( — 1 + In ( ^ 



S2 



(4.5) 



Some more complicated integrals will be obtained in the next subsection. 
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B. Partie finie based on angular integration 



The idea is to compute the partie-finie integral by performing an angular integration, 
followed by the integration over some radial variable. In a first stage, consider an integral 
that diverges at the point 1, but converges at the point 2. According to (|0|), we need to 
compute it over the domain \ Bi{s); so it is natural to change the integration variable x 
to ri = X — yi, carry on the angular integration over dQi = dQ{ni), and then, the radial 
integration over ri = |ri| varying from s to infinity, i.e. 

//• r+oo !• 

d^^F = / dhi F= dnrj I dVti F . (4.6) 
Jrx>s Js J 

In the more general case where the integral is simultaneously divergent at the two points 
1 and 2, this method stricto sensu is no longer valid since the radial integration in ( [4.6|) 
becomes divergent when ri = ri2. Yet, still it is advantageous to dispose of a mean to 
change the variable x into ri in order to obtain a convenient radial integration (even at 
the price of breaking the symmetry between the points 1 and 2). We shall derive here two 
Propositions, based on this idea, whose implementation in practical computations constitutes 
a very efficient mean to determine the partie finie, without any a priori restriction on the 



form of integrand as in the application of the Riesz formula ( [4.4| ). 

As a matter of fact, in the first proposition, the computation of a partie-finie integral 
with two singularities 1 and 2 boils down to the computation of a partie-finie integral with 
singularity 1 and a finite-part integral (FP) whose singularity is located at infinity: ri = 
|x ~ yil ~^ +00 (so to speak, the singularity 2 is "rejected" to infinity). 

Proposition 1 For any function F in the class T we can write : 

Pf,,,,, l'd'^F = Pfs^FP f d'T^ [^"j [f- J2 ^2fb] I , (4.7) 

J L ->0 J \"S2/ b+3<0 ^ ^ 

where the 2fb's denote the coefficients of the expansion of F near r2 = 0. 

In words, in order to compute the partie finie one can (i) "regularize" F around the point 2 
by subtracting out from it the terms yielding a divergence at 2, i.e. 



F2 = F - 5^ /, , (4.8) 

fe+3<0 

and (ii) compute the integral of the regularized F2 using the partie finie around 1 and the 
finite part when i? — > to deal with the divergency at infinity. Notice that the latter 
divergency has been introduced simply because of the term corresponding to 6 = —3 in ( |4.8| ) 
if non-zero. By finite part when 5 — > we mean the zeroth-order coefficient in the Laurent 
expansion of the analytic continuation with respect to the parameter 5 G C. The analytic 
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continuation is straightforwardly defined from the domain of the complex plane ^{B) > 
in which the integral converges at infinity. 

Proof. We consider two open domains Vi and V2 that are supposed to be disjoined, 
T^if]1^2 = 0; complementary in M^, i.e. T>i[JV2 = M.^, and such that yi G Vi and 
y2 ^T^2- From Definition 4, the partie-finie integral over D2 reads as (for small enough s) 

Pf / rf3xF = lim|/ £^F+ ^ [ dn2fb + \n(-) [ dn^f-s} . 

Now, two short computations reveal that 

r „b+3 r 

J2 ^2 /fe = - E / fb , (4.9a) 



b+3<0 



X \ 1 „ , I S 



Furthermore, since the integral appearing in ( |4.9a| ) is convergent at infinity, one can add 
without harm the same finite part operation when i? — > as in ( [4.9b| ). Thus, the integral 
over V2 may be re-written as 



lim If F - y FP /" d^x f A 



= lim Iff [ f ^2 - y FP [ d^^ f /i 

= lim ( FP /" c/=^x f F2 - y FP /" rf^x f /i 

We have used the facts that the integral of F converges at infinity (first equality) and 
the integral of F2 converges at the singularity 2 (second equality). Adding up the other 
contribution extending over T>i, we readily obtain the complete partie finie as 

lim I FP / d^^ f ^2 + y ^ / dn, fa + In ( -] [ dQ, 

Since the coefficients ifa, for a < —3, are those of the expansion when ri ^ of F as well 
as of F2, we recognize in the expression above the partie finie (with respect to 1 only) of the 
integral of the regularized function F2. Hence the intermediate expression 

Pfs„s2 Jd^xF = Pf.ijFP^ j ^^'x (^M^ . (4.10) 

To establish the proposition it remains to change of variable x into ri. At that point, we 
must be careful, because under this change of variable the regularization factor |x|^ changes 
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itself in a complicated way. Fortunately, we can limit ourselves to the case where B is 
infinitesimal, since we shall take the finite part afterwards, making B ^ 0. We substitute to 



x| in the right side of ( [4.1U ) its equivalent expression in terms of ri and where we expand 



when _B ^ 0, i.e. 



iB 



Bin 



rf<|l + |ln 



> "iyi 

ri 



y? 



+ 0{B' 



(4.11) 



where rii.yi denotes the usual scalar product on (and = yi.yi). Now, the dominant 
term in the latter expansion amounts simply to replacing |x|^ by rf , which would yield 
precisely the result ( [4. 71) we want to prove; but we have still to show that all the extra terms 
in the expansion (|4.11| ), which carry at least a factor B in front, do not contribute to the 
final result, i.e. that 



FP 



B 
~2 



In 



1 + 



2 "i-yi 

ri 



+ 



yf 

72" 



F2 + 0{B' 







(4.12) 



Because of the factor B in front, the only possible contribution to the finite part for 5 — * 
occurs when the integral develops a pole at i? = due to the behaviour of the integrand at 
infinity (ri — > +oo). Hence, as indicated in ( |4.12 ), the value of the integral depends only on 
the bound at infinity [this is also why we did not write a Pf^j symbol in ( [4.121) : the partie 
finie deals with the bound ri = 0, which is irrelevant to this case]. In order to evaluate the 
pole, we replace the integrand by its expansion when ri +oo. We know that F behaves as 
o^^j at a maximum |x| — +oo to ensure the convergence of the integral of F at infinity, 

swhen o we have F = o^^j when ri — > +oo. Now, from the defining expression ( [4.8| ) of 
F2, we obtain 



/-3(ni) +0 



when ri +00 , 



(4.13) 



after making the replacements of r2 and 112 by ri and ni which are permitted because 
we are working at the dominant order when ri — > +00. On the other hand, we have 

. So that the integral to be computed (as concerns the 



In 

on 



1 + + 4 

n rf 



2^^ + 

n 



y relevant bound at infinity) reads as 



d^Ti rf In 



1 + 



, yf 



Fo = -2 



+00 



dri r-^ 



B-2 



dn, ni.yi /_3(ni) + o(r5') 



This integral cannot generate a pole at i? = since such a pole could come only from a 
radial integral of the type /^°° drirf~^ (after the angular integration has been performed). 
Repeating the same reasoning to any higher orders in B, we prove the equation ( [4. 12]) as 
well as Proposition 1. 

In practice. Proposition 1 is used with the integration with respect to rii, followed by 
the integration over ri varying from (Pfsi takes care of this bound) to infinity (where 
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FPb^o does the work); Proposition 1 justifies this process even when the original integral is 
divergent at both singularities. The result of the angular integration depends on where the 
field point is located, either inside the ball Bi{ri2) centered on yi and of radius ri2 (the point 
2 lies on the surface of this ball) or in the complementary domain \ Bi{ri2). Therefore, 
a natural splitting of the integral ( [4.7D is 

\ B 

n 



Pf,,,,, d^^F= Ffs, / dhi F2 + FP / dh^ [-] F2 , (4.14) 

J JBi{ri2) -S^O iM»\Bi(ri2) V / 

taking into account the fact that the partie finie Pf,, applies only to the inner integral, over 
Bi{ri2), and the finite part FP^^^o only to the outer one, over \ ;Bi(ri2). To be more 
specific, the angular integral of F2 defines two angular- average functions l2{i^i) and J2{i^i) 
depending on whether x is in Bi{ri2) or its complementary: 



An 




dill ~ ) hin) when ri < ri2 , 

F2= ■{ _ (4.15) 
when ri > ri2 . 



The functions I2 and J2 depend also explicitely on the source points yi and y2. [As an 

^ ^ 2 - — -2 — 

example, in the case F2 = r2, we find I2 = + and J2 = ri + ^.] Now, knowing I2 
and J2, we can achieve the radial integration according to the formula 



//•r-12 ^ r+oo / \B ^ 

d^x F = Att Pf,, J dn r\ I2 + An FP j dr, i^yj r\ J2 . (4.16) 



The first term in ( |4.16| ) is quite simple to handle in practice, whereas the second one is 



more difficult because it requires a priori the knowledge of a closed-form expression for the 
integral of rf """^ J2, valid for any B such that ^{B) > 0. Obtaining this may not be feasible if 
F is too complicated; in this event, we should use a different form of the integral at infinity. 
The second proposition, which provides the appropriate form, constitutes, perhaps, the most 
powerful way to compute the partie finie in rather complicated applications. 

Proposition 2 The partie finie of the integral of F E J-" (if convergent at infinity) reads as: 
d^^ F = An Pf,, / dr, hi^r^) + An / rfn J2(ri) + - (r^F)^ 

Jo Jrri '"l 

^An{rlF)^\n(^-^ , (4.17) 
(and similarly by interchange of 1 and 2). 

Proof. Consider the angular average of the expansion of F2 when vi ^ +00 which has been 
determined in ( |4.13|) . We get 



J.^ l'^h = ~{rlF),+o(^l,] . (4.18) 
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where the coefficient of the dominant term is made of a Hadamard partie finie at point 2. 
Let us subtract and add to J2 inside the second integral in ( |4.16| ) the previous dominant 
term at infinity. In this way, we may re-write it as the sum of a convergent integral at 
infinity on one hand, to which we can then remove the finite part prescription, and a simple 
extra integral on the other hand. Namely, 



dri 



rlJ2 + -{rlF). 
ri 



{tIF\ fp 



The extra integral is finally computed in a simple way as 



dri f^iV 



FP 



±1 (HY =FF 
ri VS2/ B->o 



'B\s2J 



ln(!ii 

S2 



where we used the properties of the analytic continuation. QED. 

Thanks to Proposition 2 we are now able to compute many integrals which could not be 
deduced from the Riesz formula ( [4.41 ), unlike for ( f4.5| ). For instance. 



Pf 



Sl,S2 



Pf 



Sl,S2 



d^x Arc 
rf (ri + ra) rfg 

rfr^{ri + r2 + ri2) r 



ln(^ 

Sl 



ln(^ 

S2 



- In 2 + - 

3 3 



27r 



12 



In 



ri2 

Sl 



+ In 



ri2 

S2 



TT 



-4 



(4.19a) 
(4.19b) 



The result for the integral ([4.19b| ) is in agreement with the one that follows from a recent 
generalization of the Riesz formula to include arbitrary powers of ri + r2 + ri2, which has 
been obtained by Jaranowski and Schafer (see the appendix B.2 in |jl2|). In any case, the 
dependence of the partie-finie integral on the two constants si and S2 is given by 



+ terms independent of Si, S2 



(4.20) 



V. PARTIE FINIE OF POISSON INTEGRALS 

In this section we investigate the main properties of the partie finie of Poisson integrals 
of singular functions in the class JF. We have in view the application to the post-Newtonian 
motion of particles in general relativity, since the post-Newtonian iteration proceeds typically 
through Poisson (or Poisson-type) integrals. Consider a fixed ("spectator") point x' G 
and, for each value of x', define the function S'x'(x) = F(x)/|x— x'| where F G JF. Clearly, for 
any given x', the function S^' belongs to the class J-'\oc, introduced in Section II, Definition 2. 
In addition, when the spectator point x' coincides with the singular point yi (and similarly 
for y2), we have Sy-^ G JF. Since (as already mentionned) Definition 4 can be extended to 
functions in the class J-'\oc, we can consider the partie-finie integral 
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This is, indeed, what we shall call the "Poisson" integral of F. In particular, when the 
spectator point x' is equal to yi, we shall write 

P(y,) = Pf / d'^S.A^) = Pf / ^i^(x) . (5.2) 

The Poisson integral is not continuous at the singular point yi because -P(x'), when initially 
defined for x' 7^ yi, admits an expansion that is singular when x' tends to yi. In the 
present Section, our aim is to understand the limit relation of the integral -P(x') when 
r[ = |x' — yi| —^0, and to connect it with the "regularized" integral -P(yi) given by ( p. 21 ). 
In particular, we shall show that the "partie finie" (in an extended Hadamard's sense) of 
P(x') at x' = yi is related in a precise way to P(yi). Let us make clear straight away that 
P(x'), as a function of x' different from yi (and y2), does not belong to the class as the 
Poisson integral typically generates logarithms in the expansion when r[ 0. In particular, 
the coefficient of zeroth power of r[ in the latter expansion contains a priori a Inr^ term, 
and its partie finie in the sense of Definition 3 is in fact not finite at all, because of the 
presence of this formally infinite constant Inr'^ = —00. A possible way to deal with this 
problem, followed by Sellier in ||^, is to exclude the Inr^ (and any higher power of Inr^) 
from the definition of the partie finie. On the other hand, in applications to the physical 
problem, the constant Inr^ can be viewed as a "renormalization" constant, which is better 
to keep as it appears all the way through the calculation. Therefore, we simply include here 
the renormalization constant Inr'^ into the definition; but, for simplicity's sake, we stick to 
the name of "partie finie" in this case (although the Inr'^ makes it formally infinite). Thus, 
for a function like P admitting a logarithmic expansion: 

eN , P(x') = J2 ^r(lnr'i)^ /a,p(n'i) + o(^r[^^ when r[ , (5.3) 



a<N 
p=0,l 



we define the Hadamard partie finie of P at 1 by 

(^)i = / ^ [/o,o(n;)+ /o,i(n;)lnr;] . (5.4) 



Theorem 3 The Hadamard partie finie at 1 (in the previous sense) of the Poisson integral 
of any F E J-" reads as 

{rlF), , (5.5) 

with r[ = |x' — yi|. Furthermore the constants Si cancel each other from the two terms in 
the right side of ( \5. 3i ) (so the partie finie depends on the two constants lnr[ and lns2y'. 



fx 

ri 



P(x) + 



ln(^ 
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In words, the partie finie of the Poisson integral at 1 is equal to the regularized integral 
P(yi), obtained from the replacement x' yi inside the integrand of -P(x'), augmented by 
a term associated with the presence of the (infinite) constant Inr^. 

Proof. The fact that the constants Si cancel out (so Si is "replaced" by r[) is a trivial 
consequence of the dependence of the partie finie on si and S2 determined in (^4.201) . For our 
proof, we need the explicit expressions of the objects -P(x'), when x' is different from yi and 
y2, and P(yi), following from Definition 4. For x' 7^ yi and ri = |x — yi| ^ 0, we have the 
expansion 



= E y^^^f^) E^i^'^i f (5-6) 

(and idem 1^2), where r[ = |x' — yi|, 9£ being the multi-spatial derivative acting on x'. 
From (|3.1|), we get the expression (for x' 7^ yi and y2) 

P(x')= - — lim| / T^^F 
4vr ^-0 tiM3\Bi(s)uB2(.) |x - x'l 

l>0 a+l+3<0 ' ' -J \ 1/ J 

+ 1^2}. (5.7) 



Applying the recipe ( |5.4|) , we start by computing the angular integral over n'^^ = (x' — yi)/r^ 
(for a fixed r[) of P(x') in the form given by (|5.7| ), and consider the limit r[ — >■ afterwards. 
Since s is fated to tend to zero first, one can choose s < r[, and as we are ultimately 
interested in the limit r[ 0, we also assume r[ < ri2. To compute the angular average of 
the divergent terms in ( |5.7| ), we make use of the identities 



4:71 \ r[ J r-^ 



41^) . gj± ) (5.8b) 



4vr ^"2/ \''i2 

(where 6qi denotes the Kronecker symbol). On the other hand, the relevant formula to treat 
the integral in the right side of (|5l7|) is 



1 afn<.;), ^^^^ 

(if r[ < ri) . 



4:71 |x — x'l I _L 
n 



We split this integral into three other ones, the first of them extending over the "exterior" 
domain M.^\Bi{r[)UB2{r[), and the two remaining ones over the ring-shaped regions i3i(r'^)\ 
Bi{s) and B2{r[) \ B2{s). Hence: 
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4:71 A-n s^o 



3\Bi(ri)uB2K) n 

— / d?iLF+ I F 

'"l J Bi{r'^)\Bi{s) J B2{r[)\B2{s) ^1 



(5.10) 



Next, supposing that r[ is small enough, we may replace F in the second and third terms 
by its own expansions around 1 and 2 respectively. We find that the divergent terms in s 
cancel out, so we are allowed to apply the limit s —>■ 0. This yields 

+ ^ [ / ^^^h + (^) / /- + / 

l>0 \ i^/ b+i+3<o ' ' \ ^/ ^ 

+ oK°)} (5.11) 

(the remainder dies out when r[ — >■ 0). Under the latter form we recognize most of the terms 
composing the integral -P(yi). Indeed, we have, respectively when ri — >• and r2 — 0, 

Sy,{^) = Y,rr'faM, (5.12a) 

a 

5y.(x) = J2 nr^4—) E^2^'^2/.(n2) . (5.12b) 

Now, using the form ( p.3|) of the partie finie with the change of notation s' = r[, we find 



+ E ^-if^ii) [ E / A + (^) / „t 

+ o(ri°)}. (5.13) 

We finally evaluate the difference between ( |5.11| ) and ( |5.13| ) and look for the partie finie in 
the sense of (|5.4|) (i.e. keeping the Inr'^ term). We obtain 
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{P)i-Piyi) 



ln( ^ 

Sl 



1 



Air 1 



(QED). 



(5.14) 



The same type of result can be proved for the partie finie of the "twice-iterated" Poisson 
integral defined by 



g(x') 



Air 



Pf / £x Ix-x'lF(x) . 



(5.15) 



We find, analogously to (|5.5|) , that 



(Q)i = -— Pf / rf^x riF(x) + 



(5.16) 



For the parties finies of the gradients of the Poisson and twice-iterated Poisson integrals, we 
get 



{^^P) 



(9.g)i 



— Pf / d'^ 4 ^(^) + In f- ) {n\nF), 
An J rf \si' 



An 



Pf / rf^x n\F{^) - 



ln( ^ 



(5.17a) 
(5.17b) 



Those results are proved in the same way as in Theorem 3 (with similar cancellations of the 
constants Si). 

VI. PARTIE FINIE PSEUDO-FUNCTIONS 



A. A class of pseudo-functions 

The concept of Hadamard partie finie of the divergent integral of functions F G yields 
a natural definition of a class of pseudo-functions PfF ("partie finie" of F), namely linear 
forms on a subset of JF, of the type G G — >< PfF, G >G R, where the result of the action 
of PfF on G is denoted using a duality bracket <,>. 

Definition 5 For any function F E J-" we define the pseudo-function PfF as the linear 
functional which associates to any G E J-', such that FG = o(|x|^'^) when |x| +oo, the 
partie-finie integral of the product FG, i.e. 



< PfF, G >= Pf / rf^x FG 



(6.1) 



where the partie-finie integral is defined by ( p.l|) . 



As we can see, the pseudo- function PfF is not a linear form on JF itself but on the subset 
of JF such that the integral converges at infinity. For simplicity's sake we will always say 
that statements like ( |6.1|) are valid VG G JF, without mentioning this restriction. Note 
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also that the partie-finie integral depends on the two constants si, S2 G M"*"*, and so is the 
pseudo-function which should indeed be denoted Pfs^^^j-^. In our simplified notation we 
omit indicating si and S2- 

An evident property of the duality bracket is its "symmetry" by exchanging the roles of 
the two slots of the bracket, namely: 

V(F, G)eT\ < PfF, G >=< PfG, F > . (6.2) 

Also evident are the properties: 

< PfF, GH >=< PfG, FH >=< Pf (FG), H >=< Ff{FGH) , 1 > . 

In the following we generally do not distinguish between the two slots in <,>. Accordingly 
we define the object 

< F,PfG' >=< PfG',F > . 

Even more, we allow for a bracket in which the two slots are filled with pseudo-functions. 
Thus, we write 

< PfF, PfG >=< PfF, G >=< PfG, F > , 

which constitutes merely the definition of the new object < PfF, PfG >. 

We denote by JF' the set of pseudo-functions PfF, when F describes the class JF, intro- 
duced by Definition 5 : JF' = |pfF; F G Jf|. Later we shall extend the definition of JF' to 
include the "limits" of some pseudo-functions. Roughly, the set JF' plays a role analogous to 
the set V in distribution theory ^ , which is dual to the class T> of functions which are both 
C°°(]R^) (for what we are concerned about here) and zero outside a compact subset of M^. In 
distribution theory the set V is endowed with the Schwartz topology : a sequence (v^n^neN 
of elements of V converges to zero if and only if (i) 3nQ G N and a compact K of M.^ such 
that Vra > uq, supp(v9„) C K, and (ii) for any multi-index L = iii2 ■ ■ - iu diVn converges 
uniformly to zero. V is the set of linear forms on T> that are continuous with respect to that 
topology. In this paper we shall not attempt to define a topology on the class JF, and shall 
limit ourselves (having in view the physical application) to the definition of the algebraic 
and differential rules obeyed by the pseudo-functions of JF'. However we can state : 

Lemma 3 The pseudo-functions oj T' , when restricted to the set T> of G°°{E?) functions 
with compact support, are distributions in the sense of Schwartz : 

PfF \ e V . (6.3) 

I'D 

Proof. All we need to check is that the pseudo-function PfF i is continuous with respect to 

I'D 

the Schwartz topology 0. Consider a sequence G P tending to zero in the sense recalled 
above. Applying the partie-finie integral in the form (|3.3| ), we get (Vs' ^ 1 and VA^ G N) 
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< PfF I ,ipn> = / rf^X Fipn 

'i' i]R3\Bi(s')UB2(s') 



/-3 



Z>0 a + I+3<iV 

+ 1 ^2 + o(/'^) . 

Since and all its derivatives tend uniformly towards zero in a given compact K, 

clearly so does the sequence of real numbers < PfF i which shows that PfF i is 

Ie) I'D 

indeed continuous (QED). 

Definition 6 The product ("■") of F E T and of PfG G , and the product of two 
pseudo-functions PfF and PfG, are defined as 

F . PfG = PfF . Pf G = Pf (FG) G r . (6.4) 

In particular F . PfG = G . PfF. 

In the following, we will remove the dot indicating the product and write indifferently 

F PfG = G PfF = Pf (FG) = PfF PfG = FG Pf 1 . (6.5) 
Notice that from the symmetry of the duality bracket we have, G JF, 

< G PfF, if >= Pf / d^x FGH =< PfF, GH > . (6.6) 



Therefore, when applied to the restriction of pseudo-functions to V, the product of Definition 
6 agrees with the product of a distribution and a function ip G G°°(]R^), i.e. 

WipeV, <^/^PfF| ,(^>=<PfF| ,V<^> . (6.7) 

\t> I'D 



B. A Dirac delta-pseudo-function 

Consider, for e G M"*"*, the Riesz delta- function that we introduced in ( |3.8| ). Since 
eSi G we can associate to it the pseudo- function Pf^^i. Now, Lemma 2 [see (p^)] can be 
re-stated by means of the duality bracket as 

lim < Pf,5i,F >= (F)i . (6.8) 

This motivates the following definition. 

Definition 7 We define the pseudo-function PiSi by 

yF eJ^, < Pf5i, F >= (F)i . (6.9) 

We then extend the definition of the set T' to include this pseudo-function: Pf6i G J-'' . 
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Obviously PfSi can be viewed as the "limit" [but we have not defined a topology on JF] of 
the pseudo-functions PfeSi when e 0. The restriction of Pf5i to V is identical to the 
usual Dirac measure, 

Pf5i I =5i = 5(x-yi), (6.10) 

so that the pseudo-function Pf^i appears as a natural generalization of the Dirac measure 
in the context of Hadamard parties finies. In the following, we shall do as if 6i would belong 
to the original class of functions JF, writing for instance 

< PfF, 6i >=< Pf5i, F >= . (6.11) 

Of course, this equation constitutes in fact the definition of the bracket < PfF,6i >. 

Definition 8 For any F & the pseudo-function Ff{F6i) is defined, consistently with the 
product ( |6.4| ), by 

VGg^, <Pf{F6i),G>={FG)i. (6.12) 

We include into T' all the pseudo-functions of this type: Pf(F5i) G J-'' (that is, we consider 
•^new = ^' + + ^^2', (ind we henceforth drop the "new"). 

Notice that an immediate consequence of the "non-distributivity" of the Hadamard partie 
finie, namely {FG)i ^ (F)i(G')i, is the fact that 

Pf(F5i) ^ (F)i Pf(5i . (6.13) 

As an example, we have {r\)\ = 0; but Pf(ri5i) is not zero, since < Pf (ri^i), 1/ri >= 
1 for instance. The pseudo-function Pi{FSi) represents the product of a delta-function 
with a function that is singular on its own support, whereas this product is ill-defined 
in the standard distribution theory. However, this object, as seen as a distribution, i.e. 
when restricted to the class V of smooth functions with compact support, does exist in 
the standard theory. Using the Taylor expansion when ri — >■ of any ip & T>, that is 
TjI>o T!^Wf'9LV3(yi), we obtain 



dill ^ 



< f-i , (6.14) 



< Pf(F5i)|^, ^ >= {F^)i = E /T ^^^(yi) / 4^ 
^ i>o ' 

where if^i denotes the coefficient of l/r[ in the expansion of F when ri — > 0. Notice 
that the sum in ( |6.14|) is always finite because / < — where = clq{F) is the smallest 
exponent of vi in the expansion of F (see Definition 1). From ( |6.14 ) we derive immediately 
the "intrinsic" form of the distribution Pi{F6i)\ , that is 

ID 

Pf(F50|^ = E ^ ^^^1 f-i = E ^ (^1-1 9lS, , (6.15) 
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where Bl^i denotes the /th partial derivative of the Dirac measure (and where the sums are 
finite). We have for instance 

Pf(^) =^A^i- (6-16) 

I'D 

Note also that the distribution Pf(F5i)| can be recovered, quite naturally, from the 
Laplacian (in the ordinary distributional sense) of the bracket corresponding to the "Poisson" 
integral of Pf(F5i), i.e. formed by Pf(-F5i) acting on the function x — > l/|x — x'|. For any 
given x', this function belongs to J^\oc and we are still allowed to consider such a bracket 
(see also Section V). Thus we define 

GM--^<PW).^>^-i.(_^)^ . (6.17) 

For x' different from the singularity yi, we find, using the Taylor expansion of l/|x — x'| 
around yi. 

Clearly the function G, if considered as a function of the variable x', belongs to JF. Now, 
we see from ( |6.15[ ) that the "ordinary" Laplacian of G(x') is precisely equal to Pf(F5i)| , 
namely 

A'G'i^ = E ^ (^1^^^ = Pi{F'6[)\^ . (6.19) 

Let us point out that G has no partie finie at the point 1: {G)i = 0; so, in order to compute 
its partie finie at 1, we are not allowed to replace formally x' by yi inside the defining 
expression (|6.17|) : 



= -47r(G)i ^< Pf (F^i), - >= - = /i • (6.20) 

ri \rij^ 1 

[The function G(x') is not continuous at 1, as we can easily see from its singular expansion 

(ill)-] 

Finally let us mention how to give a sense to a pseudo-function that would be associated 
with the square of the delta-function. Ve > 0, we have sSf G JF, and hence, we can consider 
the partie-finie integral of e^^-F. In the limit e ^ we get 

lim < Fi^Sl F >= lim Pf / d^x F = 0, (6.21) 

£— »0 £^0 J 

essentially because we have a square e"^ in factor which kills any divergencies arising from 
the integral. Therefore Pf6f is (defined to be) identically zero. More generally. 
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yP eJ^ , Pi{FSl) = , (6.22) 
and we shall not hesitate to write such identities as 

< Pf5i, F6i >=< 6i, Pf (F5i) >=< Ff {F6l), 1 >= . (6.23) 

Note also that 

Fi{F6i62) = . (6.24) 

VII. DERIVATIVE OF PSEUDO-FUNCTIONS 

A. A derivative operator on T 

From now on we shall generally suppose, in order to simplify the presentation, that the 
powers of ri and r2 in the expansions of F G around the two singularities are positive 
or negative integers (g Z). Our aim is to define an appropriate partial derivative operator 
acting on the pseudo-functions of the type PfF. First of all, we know (Lemma 3) that the 
restriction of PfF to "D is a distribution in the ordinary sense, so we already have at our 
disposal the derivative operator of distribution theory [||, which is uniquely determined — 
as well as any higher-order derivatives — by the requirement: 

V(^gP, <(9,(PfF| ),(^>=-<PfF| ,diip> . (7.1) 

It is clear from viewing PfF i as an integral operator acting on ip, that ( |7.1| ) corresponds 

to a rule of "integration by part" in which the "all-integrated" (surface) term vanishes. In 
particular the "integral of a gradient" is zero. This motivates the following definition. 

Definition 9 A partial derivative operator di acting on pseudo-functions of T' is said to 
satisfy the rule of integration by parts iff 

MF.GeJ^, < a,(PfF),G >= - < 9,(PfG),F > . (7.2) 

Notice the symmetry between the two slots of the duality bracket in ( [7.21 ). As an immediate 
consequence, for a derivative operator satisfying this rule, we have 

VFgJ^, < 5i(PfF),F >= . (7.3) 

Furthermore, if we assume 5j(Pf 1) = in addition to Definition 9, then 

VFgJ^, < 9i(PfF),l >= . (7.4) 

Of course, both ( |7.3| ) and ( |7.4| ) correspond to the intuitive idea that the integral of a gradient 
(in a "distributional-extended" sense) should be zero. 
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Proposition 3 The most general derivative operator on T' satisfying the rule of integration 
by parts (|7.2|) reads 



d,{FiF) = Fi{d,F) + Bi[F] er, 



(7.5) 



where Pf{diF) represents the "ordinary" derivative, and where the "distributional" term 
Dj[F] = Hj[F] + D['^'^*[F] is the sum of the general solution of the homogeneous equation, 
i.e. a linear functional Hj[F] such that 



yF,GeT, <Ui[F],G > + <Ui[G],F >=0 
and of the particular solution defined by 

1 

k>0 



(7.6) 



(7.7) 



When applied on any G E J-", the particular solution reads 



< Br'[F], G>= dQi n\ [i /_i g^, + J2 f-2-k 9k 



+ 1^2 



fc>0 



(7.8) 



Proof. We replace the form ( |7.5| ) of the derivative operator into the rule ( [7.2[ ) and find 

< D^F], G > + < DJG], F >= - < Pf((9iF), G > - < Pf(5iG), F > . 
The right-hand side can be readily re- written as the partie-finie integral of a gradient, 

< Di[F],G> + < Bi[G],F>= -Pf j c/^x d^{FG) . 



(7.9) 



Now we know from ( p.4| ) that the integral of a gradient is equal to the partie finie of the 
surface integrals around the singularities when the surface areas shrink to zero; thus 

< Bi[F],G > + < D,[G], F >= 47i{n\rlFG)i + 1^2. 

We replace into the right side F and G by their expansions around 1, and after an easy 
calculation we arrive at 



< B,[F],G> + < D,[G],F>= [ dQ, n\ /_i + V(/-2-fc 9k+ fk 9-2 

I 11 ^-^ \i 111 

L A;>0 



1 



(7.10) 



It is clear that the particular solution given by ( |7.7| ) or ( [7. 81 ) solves the latter equation. As a 
consequence, the most general solution is simply obtained by adding the general solution of 
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the homogeneous equation, i.e. (|7.10|) with zero in the right side, which is precisely a Hj[F] 
satisfying the "anti-symmetry" property < Hj[-F], G > + < H^G], -F >= 0. QED. 

As we see from Proposition 3, the rule of integration by parts does not permit, unlike in 
the case of distribution theory [see ([7.1D], to fully specify the derivative operator. Obviously, 
we must supplement the rule by another statement indicating the cases for which the new 
derivative should reduce to the "ordinary" one, i.e. when we should have diiViF) = Pf {diF). 
Clearly, we would like to recover the ordinary derivative in the cases where the function is 
"not too much singular". In the following, we shall require essentially that our derivative 
reduces to the ordinary one when the function F is bounded near the singularities [in addition 
of belonging to C°°(R^ — {yi^2})], in the sense that there exists a neighbourhood A/" containing 
the two singularities yi and y2 and a constant M G M"^* such that x G A/" ^ 1-^(^)1 ^ 
Let us refer to the coefficients of the negative powers of ri and r2 in the expansions of F, 
i.e. the i/_i_fc's and 2/-i-a:'s where k E N, as the singular coefficients of F (recall that 
we assumed that the powers of ri and r2 are integers). Clearly, a function is bounded near 
the singularities if and only if all its singular coefficients vanish. This means that we shall 
require that the distributional term Dj[F], which is a linear functional of the coefficients in 
the expansions of F, should depend only on the singular coefficients if^i^k and ^f-i-k of 
F. This is already the case of our particular solution D^'^'^*[F] in ( |7.7| ). We now look for the 
most general possible Hj[F] depending on the i/_i_fc's (and 1 2). 

All the singular coefficients admit some spherical-harmonics or equivalently STF expan- 
sions of the type (|2.8|)-(pT9D, with STF-tensorial coefficients if^i_k [where L = ii ■ ■ ■ ii; see 
( [^.8D for definition], so we are led to requiring that Hj[F] be the most general (linear) func- 
tional of the STF tensors if^i_f. and 1^2. Moreover, we demand that HJF], like D^'^'^*[F], 
is proportional to the Dirac pseudo-function Pf(5i (as we shall see, the gradient of Pf5i is 
itself proportional to Pf5i so there is no loss of generality). Now, we have also to take into 
account the fact that the dimensionality of Hj[F] should be compatible with the one of 
Pf{diF). Endowing with a unit of length to measure the space coordinates, the Dirac 
pseudo- function Pf^i takes the dimension of the inverse cube of a length, and Hj[F] the 
dimension of F divided by this length (in physical applications, we do not want to introduce 
any special physical scale). We conclude that Hj[F] must be of the general form: 

fc>0 /=0 

where the ak/s and Pk/^ denote some purely constant numerical coefficients (and where, 
as usual, the sum over k is finite). Applying this Hj[F] on any G we readily obtain 

+ 00 

<HaF],G>=5^ J2 

k>0 1=0 

At last we must impose the anti-symmetry condition ( |7.6| ). For any G whose all singular 
coefficients vanish we have < Hj[G],-F >= 0; then, the anti-symmetry condition tells us 



ak,i nf t 



iL 

-1- 



^5i] +1^2 



(7.11) 



^+1 

21- 



k9- 



l+k 



kJ 



f iL 

7-1- 



k 9-\- 



+ 1^2. (7.12) 
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that ( |7.12|) should be identically zero for any such G and any F. Therefore, we must have 
ak^i = and Pk,i = whenever A; > 1, so we are left with only the coefficients ao^i and Po,i, 
and the condition ( |7.6| ) now implies 



+ 00 



V. 



(2«+l)!! \2l+3 



l+l 



1=0 



+ 1^2, 



which can clearly be satisfied only if (and only if) ^^«o,/ + /^o,z = 0. Thus, posing ai = ao,z, 
we have just proved: 

Lemma 4 The most general Hj[F] that vanishes for any bounded function F E and 
possesses the correct dimension depends only on (the STF-harmonics of) the singular coef- 
ficients and 2 /-I and is given by 



(7.13) 



1=0 



where the ai 's form a countable set of arbitrary numerical coefficients. 



[The angular dependence of the first term in ( [7.13| ) is expressed by means of the STF tensor 



n 



iL 



Equivalently we have 



<H,[F],G>=5^ 



ai 



(2/+3)! 



/=0 



(7.14) 



This expression is anti-symmetric in the exchange F ^ G as required. 

To sum up, we have obtained the most general derivative operator diiVlF) = P{{diF) + 
Dj[F] that satisfies the rule of integration by parts and depends only on the singular co- 
efficients of F. The distributional term DJF] is the sum of a "particular" solution fully 
specified by {fTTj ) or ( |7.8|) , and of a "homogeneous" solution given by ( [7.13|) or (|7.14|) in 
terms of an infinite set of arbitrary numerical coefficients a/ G M (and / G N). In Section 
VIII we shall see how one can reduce the arbitrariness of the definition of the derivative to 
only one single coefficient i^' G M. 



B. Some properties of the derivative 

At this stage, one can already investigate some properties of the distributional term 
I)i[F] = Df'*[F] + Hi[F], using the fact that the yet un-specified < Hi[F],G > depends 
only on and ig^i (and 1^2). Let us first check that the derivative operator, when 
restricted to the smooth and compact-support functions of V, reduces to the distributional 
derivative of distribution theory [Q. This must actually be true since the fundamental 
property ( [7.1] ) of the distributional derivative is a particular case of our rule of integration 
by parts, and because the derivative of G reduces to the ordinary one. However, it is 
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instructive to verify directly this fact using the expression ( [7. 7|) . Applying Dj[F] on G P 
and using the Taylor expansion of ^ around 1: (y9 = X]fc>o we obtain 

Hence the intrinsic expression of the distributional terms on T>, 

DJ^] I = E ^^^^1 / f-2-k + 1^2, (7.15) 



A;>0 



which agrees with the distributional part of the derivative of a function with tempered 
singularities in distribution theory. For example, we can write 



(7.16) 



However, when acting on functions of the full set JF, the derivative generally leads to 
properties which have no equivalent in distributional theory. For instance, although the 
distributional derivative of l/rf reduces on V to the ordinary derivative, i.e. Di[^] I ~ 0> 
on JF it does not: 



dAPi\] =Fi( -2^^+471 n\5i] . (7.17) 



For the distributional derivative of 1/r^ on JF we find 

a.(Pf;^)=Pf(-3^ + 4-^*0' 

The expression of the distributional term is apparently different from the corresponding 
result ( |7.16| ) in distribution theory. However we shall see after learning how to differentiate 
the Dirac pseudo-function Pi6i that the distributional term DJ^] takes in fact the same 
form on as on X> [see ( 7.28| ) below]. 



We come now to an important point. In this paper we have defined a "pointwise" product 
of pseudo-functions (see Definition 6), which reduces to the ordinary product in all the cases 
where the functions are regular enough. For instance, it coincides with the ordinary product 
for C°° functions, or even continuous or locally integrable functions (adopting the class J-'ioc) ■ 
Next, we introduced a derivative operator that acts merely as the ordinary derivative for 
a large class of not-too-singular functions (those which are bounded near the singularities, 
see Proposition 3). In particular, the derivative is equal to the ordinary one when the 
functions are at the location of the two singularities. However, we know from a theorem 



of Schwartz that it is impossible to define a multiplication for distributions having the 
previous properties and such that the distributional derivation satisfies the standard formula 
for the derivation of a product (Leibniz's rule). In agreement with that theorem, we find 
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that the derivative operator defined by (|7.5|) -(fr7| ) does not obey in general the Leibniz rule, 
whereas it does satisfy it by definition in an "integrated sense" , namely: 



< d,[Pf{FG)], 1 >= =< di{PfF)G + Fdi{PfG), 1 > 



(7.19) 



However it does not satisfy the Leibniz rule in a "local sense" , i.e. we have, generically for 
two functions F, G G JF, 



a,[Pf(FG)] - <9i(PfF)G - F9,(PfG) . 
This means that, a priori, 

< di[Fi{FG)], H> - < a,(PfF), GH> - < (9i(PfG), FH>^0, 
or, equivalently, since the Leibniz rule is satisfied by the ordinary derivative, 
< Bi[FG],H > - < Bi[F],GH > - < Bi[G],FH . 



(7.20) 



(7.21) 



(7.22) 



Actually, in accordance with the theorem in , ( |7.20| ) must be true even when the pseudo- 
function is regarded as a distribution on V. To check this, let us compute the left side of 
( [7.22| ) in the case where Dj is the particular solution D^*^"^* defined by ([7.7[), and where H is 
equal to some cp E T>. We employ the Taylor expansion of around 1 and 2, and, strictly 
following the definition of the distributional term in (|7.7| ), we arrive at 

= V ^OkS, [ dn, n\nf h /_i + i - V g,^k-i +1^2. (7.23) 

^-^ kl I 1 1 ^11 1 1 J 



k>l 



j=0 



The right side of ( |7.23|) equals ^di6i in the case where F = ^ and G = for instance. It is 
not possible to add a homogeneous solution of the form ( [7.13|) so as to get always zero. As 
the result ( [7.23| ) depends only on the singular coefficients of F and G, we recover the Leibniz 
rule whenever F or G is bounded near the singularities. Besides, we can verify directly on 
( [7.23| ) that the Leibniz rule is indeed true in an integrated sense, since the integral over 
of ( [7.23[ ) picks up only the term with k = which gives no contribution. 



C. Derivative of the Dirac pseudo- function 

In this subsection we compute the distributional term < Dj[F],G > given by the sum 
of ( |7.8| ) and (|7. 14|) assuming that either F or G is equal to the Riesz delta-function = 
^^l^rf"'^ for some small e > 0. (We come back for a moment to Definition 1 in which the 
powers of ri and r2 in the expansions of F or G are real.) We notice first that the terms 
depending on the singular coefficients and ig_i are present only when the exponent — 1 
belongs to both families of indices (aj)^^^ corresponding to F and G (remind Definition 1). 
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This means that, choosing e to be different from 2, these terms will not contribute to the 
present calculation, and in particular that the homogeneous part < Hj[F],G' > will always 
give zero, provided that either F or G is equal to ^Si. From the expression ( [7.8| ) we get 

< D,,[,5i], G> = 6(1-6) J^n\ g,_, , (7.24a) 

< DJF], ,5, > = e(l - ^) ^ d^r{f I nf f^^-i • (7.24b) 

Furthermore, by choosing 6 smaller than the spacing between some exponents at of G 
(specifically e < 1 — with a^^ is such that a^^ < 1 < aj^+i) we can arrange for having 
iQi^e = so that ( [7.24a| ) becomes identically zero. Anyway, in the limit e — we come up 
formally with both relations < Dj[5i],G' >= and < Dj[F],5i >= 0. The former tells us 
that the distributional derivative of Pf5i reduces to the ordinary one, i.e. 

9,(Pf5i) = Pf(5,5i) . (7.25) 



The latter [that we already knew from ( |6.23| )] shows via the rule of integration by parts that 



the action of di(Pf6i) over any function F E is equal to minus the action of Pf^i over the 
derivative diF. 

Definition 10 The derivative of the Dirac pseudo-function Pf5i is defined by 

VF e ^ , < 9i(Pf5i), F >= - < Pf5i, diF >= -idiF)i . (7.26) 

We can summarize the properties of the derivative of the Dirac pseudo-function by writing 
the successive identities 

< a,(Pf(5i),F >=< Ff{dA),F>= - < Pf(5i,9,F >= -(9,F)i , 
as well as similar identities obtained by exchanging the roles of F and 6i, 

< 9,(PfF),5i >=< Pf(9,F),5i >= - < FfF,dA >= (5.F)i . 

Lemma 5 The intrinsic form of the derivative of the Dirac pseudo-function is 

a,(Pf5i) = -Pf(3^5i) . (7.27) 



The proof is evident from using the identity ( p.ll|) . The form ( |7.271 ) [with ( |7.25| )] is quite 



useful in practice; for instance, it permits us to re- write the derivative of the pseudo-function 
Pf(^) as computed in ( [7.18| ) into the form 
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(7.28) 



where the distributional term takes the same form as in the distribution theory [compare 
with ([7l6|) ]. 

The preceding definition and lemma are easily extended to the case of the pseudo- 
functions Pf(-F5i). The derivative of these objects is defined by the mean of the relation 



< d,[Pi{FSi)],G>= - < Pf(F5i),9,G >= -(F9,G)i . 
Then, from the identity (|2.11| ), we readily get the intrinsic form 

Notice the interesting particular case 















)] = 


0, 



(7.29) 



(7.30) 



(7.31) 



which is also an immediate consequence of ( [2.161 ). Finally, let us mention that the Leibniz 
rule happens to hold in the special case where one of the pseudo-functions is of the type 
Pf(G(5i), i.e. 



di[PfF. Pf(G5i)] = diiPfF). Pf(G5i) + PfF. di[Pf{G6i)] 
(the verification is straightforward). 



(7.32) 



VIII. MULTIPLE DERIVATIVES 
A. General construction 

From Proposition 3 we can give a meaning to 



< ^^(PfF), G >= Pf j d^x diFG+ <Di[F],G> , (8.1) 

which will be also denoted < 9j(PfF),PfG >. We now define the more complicated object 
< (9j(PfF), 9j(PfG) >. Since the distributional term DJF] has the form Pf(if5i) plus 
1 ^ 2, and because §^-{^^2^ entail such identities as < Pf (G(5i), Pf (i/5i) >= =< 



Pf (G5i), Pf (if(52) >, we deduce that the duality bracket applied on any two distributional 
terms is always zero: 

VF,GgJ^, <Bi[F],Dj[G]>=0 . (8.2) 

When constructing the bracket < di{PfF),dj(PfG) > we shall meet a product of two dis- 
tributional terms which gives zero by ( |8.2| ), and we shall be left only with the ordinary part 
as well as the two cross terms involving one distributional term. Therefore, 
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< di{PfF),dj{PiG) > = Pf y c/^x d.FdjG 

+ < Bi[F],djG > + < Dj[G], diF > . (8.3) 
[The ordinary part can equivalently be written as 

Pf J d^x diFdjG =< Fi{diF),Fi{djG) >=< Fi{d^F),djG >=< diF,Fi{djG) > .] 

We now intend to introduce the second-order derivative operator. The generahzation to any 
/th-order derivative is straightforward and will be stated without proof. By extending the 
rule of integration by parts presented in Definition 9, we are led, quite naturally, to require 
that 

VF, G G J^, < dij{FiF), G>=- < dj{FfF), di{FfG) > , (8.4) 

where the object < dj{F{F),di{FfG) > has just been given in (|8.3|) . For the moment, we 
are careful at distinguishing the order of the indices i and j. Let us look for the expression 
of the distributional term Djj[F] corresponding to the double derivative, viz 

a,,(PfF) = Pf(9,,F) + D,,[F] , (8.5) 



in terms of the single-derivative term Dj[F]. Inserting ( ^.5|) into the required property ( ^.4| 
we arrive immediately at 

< I)ij[F],G>= -Pf j rf^x di{djFG)- < Dj[F],9iG > - < D,[G],6»jF > . 



Next recall the formula (|7.9| ) which tells us that any partie-finie integral of a gradient is the 
sum of two distributional contributions. Using this property we obtain the simple result 

< D,,[F],G > = < Da9,F],G > - < D,[F],a,G > 

= < T)i[djF],G > + < diT)j[FlG > . (8.6) 



The formula ( [7.29| ) allowed us to obtain the second equality; so the intrinsic form of the 



second-order distributional term is obtained as 

D,,[F] = D,[a,F] + 9,D,[F]. (8.7) 
This result is easily extendible to any multiple derivatives, demanding that, to any order /, 
< 9,,,,...,,(PfF),G >= - < 9,,...,,(PfF),9,,(PfG) > , (8.8) 



where the right side is obtained in a way similar to ( p.3|) . We can even impose the more 
general rule of integration by parts, that for any k = 1, ■ ■ ■ ,1, 

< 9,,,,...,,(PfF),G >= < 9,,^,,,^,...,,(PfF),5,,,,_,..,,(PfG) > . (8.9) 

Then the following is proved by induction over /. 
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Proposition 4 If a multi-derivative operator 

di,i,...i,iFfF) = Fi{di,^,...i,F) + D,,,,...,JF] , (8.10) 

satisfies the rule of integration by parts ( |8.8| ) or ( |8.9| ), then the l-th order distributional term 
Dj^j2...jjF] is given in terms of the first- order Di^[F]'s by 

I 

^n^2■■■^^[F] = J] 9,,...,,_,D, J9,,^,..,,F] . (8.11) 
k=l 



Recall that this result is valid for any distributional derivative of the form given by Propo- 
sition 3, i.e. Dj[F] = D^^'^*[F] + HJF]. Therefore, the rule of integration by parts has 
permitted us to construct uniquely all higher-order derivatives from a given choice of first- 
order derivative Dj[F], i.e. from a given choice of "homogeneous" solution Hj[F]. Notice 
that a priori this construction does not yield some commuting multi-derivatives (i.e. the 



Schwarz lemma is not valid in general), because evidently the right side of the formula (^.11|) 



is not necessarily symmetric in all its indices. However, as a central result of this paper, 
we shall show now that it is possible to find an initial HJF] such that the derivatives do 
commute to any order. 

Theorem 4 The most general derivative operator di{PfF) = Pi{diF) + ^i[F] such that 
(i) the distributional term Dj[F] depends only on the singular coefficients of F, 
(a) all multi-derivatives satisfy the rule of integration by parts, 
(Hi) all multi-derivatives commute (i.e. the Dj^j2...jjF] 's are symmetric in 11X2 ■ ■ - ii), 

is given by 



+00 / iL \ 

B,[F] = 47r J] Pf ( G [nf f_\ - n\ /_\^] n^i + J] \ iVfe^i +1-2 



(8.12) 



where the coefficients C/ = (/ + 1) K + X]j=i J+T 



depend on an arbitrary constant K . 



(Actually the theorem states that the derivative operator depends a priori on two different 
constants K\ and K2 for each of the two singularities. In the following we shall assume 
for simplicity that the constants are the same, so that the way to differentiate does not 
distinguish between the different singularities.) Notice that Dj[F] differs from the partic- 
ular solution D^'^'^*[F] given by (|7.7|) only in the terms depending on the "least singular" 
coefficients and 2/-i- 

Proof. According to the assumptions (i) and (ii) we already know (see Proposition 3 and 
Lemma 4) that the distributional term must be of the form Dj[F] = Df^'^*[F] + Hj[F], where 
the particular solution is given explicitly by ( [7.7| ), and where the homogeneous term takes 
the form ( |7.13| ) depending on a set of arbitrary coefficients «/. Furthermore, we know from 



38 



Proposition 4 that all higher-order derivatives are generated from the first-order one in the 



way specified by It only remains to show that the coefficients ai can be computed 

in order that the assumption (iii) of commutation of derivatives be fulfilled, and that the 
derivative is given by ( ^■121 ). 

What we want then is to impose the symmetry of Dij[-F] in ij. We compute the anti- 
symmetric projection [ij] = of the second-order distributional term associated with the 



particular solution ( [7.7|) 



Dj:;f[F] = Df;-[9,]F]+9[,D;p[F] 



.part 



part r 



(8.13) 



The first term is readily obtained from ( |2.12| ) which tells us that the ath coefficient in the 
ri-expansion of the gradient is ifa{djF) = (a + l)n{ifa+i + i/a+i- On the other hand, 
the second term in ( ^.131 ) comes directly from using the formula ( [7.3U| ). It follows that the 
anti-symmetric projection depends only on the expansion coefficients i/o, and 1^2 
through the simple formula. 



Dg[F] = 27rPf(nf [n^/f /o+^/f A 
or, using the relation ( p.l3| ) for the operator dl, 



(5i + 1 ^ 2 



(8.14) 



Dg[F] = 2vr^(/ + l)Pf 



1=0 



5i] +1^2. 



Note that by applying this on any G, we get 



(8.15) 



<D-f[F],G>=-2.5^ 



(/+!) {l+iy. ( fL[i ;^j]L. fL[i ^j]L 



fo ' rr+ a or + 1 - 2 . 



1=0 



(2^+3)!: J ^ J ^ ^- 



Next, we add the homogeneous solution. By performing a computation similar as the pre- 
vious one (but a bit more involved) we find, based on the expression ( |7.13| ), 



^mI^] = E ^[(^ + - + Pf (^f h /o''^+ f-l%) + 1 2 . (8.16) 



1=0 



Remarkably, H[jj][F] takes exactly the same form as ( ^.15| ). Hence, we are able to determine 
a relation to be satisfied by the looked-for coefficients for any / in order that the non- 
commuting part ( p.l5[ ) associated to the particular solution be cancelled out by that of the 
homogeneous one: V/, {I + 2)ai — {I + = — 27r(2/ + 3). Given any initial value for ao 

the solution reads as 



ai = il + I] 



+ 27ry - + - — - 



-27r + 47r(/ + 1) + V 



1 



(8.17) 
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in which we have introduced the new arbitrary constant K = ^ + \- Inserting ( p.l7| ) back 
into the expression for Dj[F] leads to the announced result ( |8.12| ). At last, we find that for 
any choice of the constant K the second- derivative operator commutes, i.e. 



D[,,][F] = H[,,][F] + D|:^f[F] = 0. 



(8.18) 



Let us verify from (|8.18|) that all higher-order multi-derivative operators commute as well, 
i.e. Dj^j2...jjF] given by the formula (|8.11|) is symmetric in all its indices. This is easily 
proved by induction over /. Suppose that to the (/ — l)th order Djjj2...j;_ JF] is symmetric, 
and re-write the formula (|8.11|) into both forms 

= D,,..,.J9,,F] + 9,,..,,_,D,JF]. 

Clearly, Dj^...jjF] is symmetric with respect to both ii ■ ■ - ii-i and Z2 ■ ■ - iu so must be sym- 
metric in all its indices (the symmetry with respect to the first and last indices being a 
consequence of the other symmetries). QED. 

We should mention that the dependence upon the arbitrary constant K of the derivative 
operator defined by Theorem 4 is 



D,[F] I =AtiK y(/ + l)Pf 



ri5i -f 1 



(8.19) 



which can also be cast into the more interesting form 



Pf(r?/-i5i; 



1^2. 



(8.20) 



We see that the "ambiguity" linked with the constant K when deriving the pseudo-function 
PfF is related to an ambiguity resulting from the addition of the term — 47rfC Pf (rj i/_i5i) + 
1 ^ 2 to PfF. In a sense, one can also view the constant K as a measure of how much the 
distributional derivative of the pseudo-function Pf(l/ri) differs from the ordinary one, i.e. 



D, 



1 

ri 



AnKFi {n\ri6i) . 



(8.21) 



Indeed, for functions which are more singular than a simple 1/ri, there is no dependence on 
the constant K; see e.g. ( [7.17D -( [7rT^ ). 



B. The Laplacian operator 

Let us compute the second-derivative of Pf(l/ri) using the formula Djj[l/ri] = 
Di[—n-[/rf\ + 9jDj[l/ri]. The first term is obtained directly from the definition ( ^.12 ) 



and the second term is computed with the help of the formula ( [7. 30] ) applied on ( ^.211 ). As 
a result, we get 
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D 



1 

ri 



An 



(8.22) 



where n*/ 



i„j 1 



l"-! 3 



Evidently (because of the trace- free re*/), when we restrict ourselves 



to smooth functions of the set V, we recover the usual formula of distributional theory, 



1 

ri 



6'^ 6i . 



(8.23) 



Since the dependence over K in ( p.22| ) drops out when taking the trace over the indices ij, 
we have Djj[l/ri] = — 47rPf5i (even on the set JF). This means that the Laplacian of 1/ri 
on JF takes the same form as the well-known formula of distribution theory: 



A ( Pf- 



-Att Pf (5i . 



We infer from the rule of integration by parts that 

< A(PfF), - >=< A ( Pf- 1 , F >= -Att (F)i 



(8.24) 



(8.25) 



which can be phrased by saying that the Poisson integral of the Laplacian of a singular 
function, as evaluated at a singular point, is equal to the partie finie of the function at that 
point. More generally, the Laplacian acting on any pseudo-function in JF' is defined by 



A(PfF) =Pf(AF) + D,,[F] , 
where the distributional term is given by 

B,,[F] = d,B,[F] + B,[d,F]. 



(8.26) 



(8.27) 



Proposition 5 Under the hypothesis of Theorem 4 the distributional term associated with 
the Laplacian operator reads 



+ 00 



D,,[F]=47r^Pf (/ + l)a 



1=0 



L\ + ri /o^ 5i - Y.^2k + 1)^ f_\_,5, +1^2. 
1 1 J ^-^ rf 1 / 

fc>0 1 ^ 

(8.28) 



The proof is straightforward and will not be detailed. Note that the dependence on K occurs 
only for functions owing some non-zero coefficients or i/o, or 1 2; for instance 

D,,K] =87rirPf (nln^i) , 
D..K/ri] =87r(/^-i)Pf (ni^i) . 

But, for more singular functions like we have 



rl rf 



(8.29) 
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Lemma 6 The Laplacian of the pseudo-function Pf(F(5i) is given by 



A 



Pf(F(5i 



Pf rlA 



F 



61 



(8.30) 



The proof is similar to the one of the formula ( [7. 301) . Two immediate particular applications 
are 

6 



^(Pf^i)=Pf(;:2^i 



A 



0, 



(8.31a) 
(8.31b) 



which can also be deduced respectively from ( p.l9|) and (|2.17|) . [( ^.31a| ) is in agreement with 
( |6.16| ).] Let us add that 



A 



Pf(t. 



Pf('^ 

ri 



(8.32) 



In practice, Lemma 6 may be used to determine some solutions of Poisson equations "in the 
sense of distributions" on JF. For instance combining ( |8.31a| ) with the formula ( p.29| ), we 
can write 



A 



Pf 



1 



IOtt 



5i 



Pf 



6 



(8.33) 



which provides a solution of the Poisson equation with source Pf (6/r^) in the sense of these 
distributions. Such a solution is by no means unique, since, from Lemma 6, one can add to 
it any "homogeneous" solution of the form Pf (if'^°™(5i) where is the product of rf with 

an arbitrary solution of the Laplace equation. Notice that ( p.33|) as it stands is well-defined 
in distribution theory and so takes the same form when restricted to V {A61 is meaningful 
on this set). However, 



Pf(^^ + 67rri5i 



(8.34) 



has no equivalent in distribution theory. 



IX. TIME DERIVATIVE AND PARTIAL DERIVATIVES 

The functions F E depend on the field point x and on the two singular source points 
yi and y2. We shall now consider the situation where the two source points represent the 
trajectories of actual particles, and therefore depend on time t. We assume that the two 
trajectories yi(t) and y2(t) are smooth, that is yi,y2 G C°°(M). In general (e.g. in the 
application to the problem of motion of point-particles) the function F will also depend on 
time through the two velocities vi(t) = dyi{t)/dt and V2(t) = dy2(t)/dt. We suppose that 
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F is a smooth functional of vi and V2. Therefore, in this section, F is supposed to take the 
form 

F = F(x,t;yi(t),y2(t)) G ^. (9.1) 

We want to investigate the partial derivatives (in a distributional sense) of the pseudo- 
function PfF with respect to the source points yi and y2, as well as the derivative of PfF 
with respect to time t. Obviously, the partial derivatives idi = d/dyi and 1 <-> 2 are closely 
related to the time derivative dt = d/dt on account of the fact that 

dtF = F + vl diF + vi diF , (9.2) 

(in the ordinary sense), where F denotes the contribution of the time-derivative due to the 
dependence over the velocities, i.e. F = a\dF/dv{ + a\dF / dv\ {'A ^^^1 denoting the two 
accelerations). In applications it is frequent that F depends on the trajectories only through 
the two distances to the field point ri = x — yi and r2 = x — y2; in that case 

diF+ diF+ diF = 0. (9.3) 

1 2 

The general function does not necessarily satisfy the latter identity. However, let us 
guess from ( |9.3| ) the result for the distributional terms iDj[F] (and 1^2) associated with 
the partial derivative idi acting on the pseudo-function PfF. Since we have supposed that 
the dependence of F on the velocities is smooth, the distributional terms will depend only 
on that part of the function which becomes singular when ri — > 0, and so, because as far 
as the singular part is concerned, the function behaves like ( |9.3|) , the distributional terms 
iDj[F] and 2D4-^] should satisfy 

Bi[F]+ I)i[F]+ Bi[F] = . (9.4) 

Now, from Theorem 4, we know that Dj[F] can be naturally split into two parts associated 
respectively with the singularities 1 and 2. Therefore, we expect that the correct distribu- 
tional term iDj[F] is equal to minus that part of Di[F] which corresponds to 1. Namely, 
using (111), 

+°° / iL \ 

(and idem for 2). This expectation is confirmed by the following definition and proposition. 

Definition 11 The partial derivative ^di (and 1^2) acting on pseudo-functions is said to 
satisfy the rule of integration by parts iff 

yF,GeJ^, <di(PiF),G> + <di(PiG),F >=di <PfF,G> . (9.6) 
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Similarly, the time derivative dt is said to satisfy the rule of integration by parts iff 



Notice that < PfF, G >= Pf / d^^FG is a function of the source points yi(t) and y2(t), as 
well as t independently if either F or G depends on the velocities. The time derivative in 



the variable t occuring through yi(t) and yi(t), and the independent t coming from the 
velocities. Let us now state a result analogous to Theorem 4, whose proof will not be given 
since it represents a simple adaptation of the one of that theorem. 

Proposition 6 Under the hypothesis of Theorem 4 the partial derivative with respect to yi 
(and idem with 1 ^ 2) is determined as 



< 9t(PfF), G > + < dt{PfG), F>= — < PfF, G > 



(9.7) 




means the total time derivative we get by taking into account both 



^i{PiF)=F^{^^F)+I)^[F], 



(9.8) 



where iDj[F] is given by ( |9.5|) . And the time derivative is determined as 



dt{PfF) = Pf{dtF) + Dt[F], 



(9.9) 



where 'Dt[F] is given by 



Bt[F]=v\ B^[F]+vl Di[F]. 

1 2 



(9.10) 



Higher-order derivatives are constructed as in Section VIII. We find for instance 




(9.11) 



Idem for the second-order time derivative, which reads as 



d^^FiF) = Fm^F) + Bt[dtF] + dtBtlF] 



(9.12) 



where dfF is given by ( p.2[ ) and Dt[F] is defined in ( p.lO| ). Furthermore the mixing up of 
derivatives of different type is allowed, and proceeds in the expected way. E.g. 



didjiPiF) = Fi{didjF)+ Bi[djF]+ diBAF] • 

12 12 12 12 



(9.13) 



Another example is 



dAj{PiF) = Fmd,,F) + Bt[d.,,F] + dtB,[d,F] + dtd.,B,[F] . 



(9.14) 
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APPENDIX A: PROOF OF THEOREM 2 



Basically the proof establishes the legitimacy of commuting some discrete series with 
integrals. Consider F ^ T . We start by evaluating the integrals 




(Ala) 



and Pf.^o / V f ^) f^) rUa(ni) , (Alb) 

where the i/a's are the coefficients of the expansion of F when r\ 0, and where Si(s) is 
the ball centred on yi and of radius s G M^* (chosen to be s < ry2). From the definition 
of the class T the sums over a in ( [Al| ) are finite. When the real part of /3 is such that 
< 5?(/9) < 1, the integrand of ( |Ala) ) is majored by 



^ ' a<— 3 ^ ' ^ ' a<— 3 



Si 



which can be integrated on i3i(s). Thus the theorem of dominated convergence of an integral 
can be applied, with the result that 



a+3<0 ' -J \ i-/ J 



The second integral is more difficult to compute because the limit a — > does not commute 
with the integration sign. We must expand as a power series of ri. Vri < ri2, 



2 



rf,(l + 2n,niA + 4)''' = rf, £ Cr^/^(-n,n,,) f^V , (A2) 



where ni2 = (yi — y2)/''"i2, and where C^{t) denotes the Gegenbauer polynomial, which is by 
definition the coefficient of in the power-series expansion of the function (1 — 2tx + x^)~^ 
when X — * (with A, t G C). See e.g. Morse and Feshbach [^] p. 602. When t G M and is 
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such that |t| < 1 (as is the case here since t = — ni.ni2), we can obtain a majoration of the 
Gegenbauer polynomial. From the formula (c/ Gradshteyn and Ryzhik p. 1030) 



C,^(cosi 



>^ T{X + k)T{X + h) 



k,h>0 
k + h=l 



k\h\ [T{\)Y 



we find that V/ 7^ 0, \Cf{co?,9)\ is always less than 



E 

fc,h>0 



(|A| + fc-l)(|A| + fc-2) 



A| + /i-l)(|A| + /i-2)---|A| 



Cf'(l). 



k\h\ 



Therefore the series |C; '^^^(— ni.ni2))(ri/ri2)'| is bounded by (1 — 2r\jrYi + ^1/^12)'^'^^, 
and thus admits a limit. Thus ( |A^ ) converges absolutly and (when 3?(q;) is large enough) 
the signs / and ^ can be interchanged: 



P/2 



where ^^"^ = Ci ^^^(—111.1112). We obtain the two terms 



ri2 

S2 



p +00 



E E 



/ I , , N I d^^^ fa Ci 

1=0 .t^o s'ir[,{a + a + l + 3)J I 



-P/2 



rvA \- - ( — 

S2J ^ a\sij 



1^0 
finite sum 



-13/2 



The finite part when a ^ of the second term reads simply 

/3 / \ +00 



/3/2 



(A3) 



On the other hand, in order to treat the first term, we must justify the commutation of the 
finite part with the infinite sum. Consider the series 



+00 

E 

1=0 



a + a + I + 3 \ri2 



For a in a disk of the complex plane centered on and of radius e (with < e < 1), we can 
bound the generic term of that series (for large enough I) by 



1 



la + / + 31 - e V ^"12 



rffii fa C;"'^ 
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which is independent of a, and whose corresponding series in / converges. Therefore we can 
apply the hmit a — > through the summation over / and deduce 

hm|(^) y y — - fdo. Lc-^i^ 



^ a + / + 3 V ^12 

a+3<0 i=0 ^ 

a + i + 3#0 



Next we apply the finite part Pf/3^0 to the sum of ( |A3D and ( [A4D , which involves finding 
the limit when /5 — > of the series 



a + / + 3 V ^12 

£i+i+3#0 



In any case the absolute value of the quantity under the sign ^ is smaller than 

Furthermore we know that Cl{\) = T{2\ + /)/[/! r(2A)]. For / ^ 0, Cj^'(l) = (2|A| + 
/ — 1)(2|A| + / — 2) • • ■ (2|A|) / /! is manifestly an increasing function of |A|, and, for / = 0, 
Cg'^'(l) = 1 is constant. From this we infer that V/ and for f3 in the disk centred on and of 
radius e, < Ci^'^{l) holds, which leads to the /3-independent bound 



which is manifestly the general term of a convergent series. Therefore the series (^^) pos- 
sesses a limit when (3^0 which is simply obtained by setting (3 = under the sign 
Using Ci{cos6) = 610 we find this limit to be if a = —3 and 

dn^ fa (A6) 



a + 3 



if a 7^ —3. Gathering the results (|X3|)- (|A4|) and (|A6|), we arrive to 

o+3<0 \ i/ 

= FP.-.o/ ('jX ('-l)\ua . (A7) 



a+3<0 
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from which we can now easily prove the equivalence with the Hadamard partie finie. Like 
in the proof of Proposition 1 we consider two open domains Vi and V2, disjoined and 
complementary in M.^, and such that Bi{s) C Vi and -62(5) C 1^2- We can write 



(A8) 



where each of the objects is defined by complex analytic continuation in a neighbourhood of 
a = (3 = [proof similar to the one after (O)]. Like in ( [4.8| ) we associate to F the function 



Fi representing its "regularization" around the point 1, 



a+3<0 



and we re-write the right side of (|A8|) as 



/ #x(^)7:i)V+/ #x(^)Y^yF.+/ A(ii)Y:i)'5].;/„. 



Of these three terms, the first two are well-defined when a and (3 tend to zero, hence their 
finite parts are simply obtained by posing a = = /5. On the other hand we have proved 
previously that the finite parts Pfa-.o and Pf/3^0 of the third term are equal and we have 
found their common value to be given by ( |A7| ). This shows immediately that 



FP.-.0 / d!^^ {—] f— ^ F = I d^:siF+ I d^xFi 



a+3<0 ' >^ \ 1/ J 



(and idem with Pfp^o). We recognize in the right side of ( |A10| ) the Hadamard partie finie 

a— *0 

of the integral. Indeed the second term clearly admits an expansion in positive powers of s, 



ViV G N , / d^^F,= V ^ / dn, fa + o(s^) , (All) 



o«+3 

'^i-(^) 0<a+3<N 

SO we recover exactly the partie-finie integral over T>i in the form given by ( p.3|) . QED. 
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